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I. INTRODUCTION 



What is the scattering matrix of a matter-coupled gauge theory? In general, this is a 
hard question involving both conceptual and technical subtleties. Nevertheless, scattering 
amplitudes enjoy a much simpler structure than implied by their expansion in terms of 
Feynman diagrams. For some theories, additional off-shell and on-shell symmetries simplify 
the amplitudes enormously. The further simplification exhibited in the planar (or infinite- 
color) limit may even allow a complete answer to the question. 

The M = 4 or maximally supersymmetric Yang-Mills theory (MSYM) may be such a 
theory. The simplifications inherent in the larger symmetry have already allowed explicit 
calculations of scattering amplitudes well beyond those for other theories. At weak cou- 
pling, advances in multi-loop and multi-leg calculations jl-7| have opened the possibility of 
probing the structure of the scattering matrix to high order in perturbation theory. The 



BDS conjecture 



8| for the all-loop resummation of maximally helicity-violating (MHV) am- 



plitudes (based on an earlier relation jo| linking one- and two-loop amplitudes) provides an 
example of possible structures that can emerge. At strong coupling, the leading expansion 
of scattering amplitudes has been computed using the AdS/CFT correspondence [10| by Al- 



day, Gaiotto, Maldacena, Sever, and Vieira jl lMl4f| . This two-sided approach, together with 



the recent developments in the evaluation of scattering amplitudes at strong coupling for 
any number of external legs and the realization that the relation between certain scattering 

n 

amplitudes and null polygon al Wilson loops carries over from strong coupling [15[ to the 



weak-coupling regime 16|-|20) offer circumstantial evidence that the M = 4 super- Yang-Mills 



theory may ultimately be solvable in the planar limit. 

Gluon MHV amplitudes, with two external legs of negative helicity and the rest of positive 
helicity, are the simplest amplitudes in a gauge theory. They are particularly simple in the 
planar limit of MSYM, where they are determined by a single helicity structure, accompanied 
by a function of scalar and pseudo-scalar momentum invariants, to all orders in perturbation 



theory. This simplicity is of course shared by their parity conjugates, the MHV amplitudes. 
The structure of the remaining non-MHV amplitudes is more complicated. At one loop, 
explicit expressions are known 2l|, |22( at arbitrary multiplicity for next-to-MHV (NMHV) 



gluon amplitudes, with three gluons of negative helicity from the rest. These are already 
more intricate, with the number of independent helicity structures growing cubicly with the 
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number of external legs, each multiplied by an independent function of scalar and pseudo- 
scalar invariants. How does this structure generalize to higher loops? No explicit expressions 
are known to date for higher-loop non-MHV amplitudes. In this paper, we take a first step 
towards filling this gap, computing the parity-even part of the two-loop six-gluon NMHV 
amplitude. These results were first reported in ref. 23] . This amplitude, which comes with 
three inequivalent helicity configurations, is the simplest non-MHV amplitude. 

General results on the structure of infrared divergences in massless gauge theories suggest 
on one hand, that the divergent terms have a simple iterative structure; and on the other, 
that all planar amplitudes with fixed number of external legs share the same structure 
of infrared-divergent terms. Together with the structure of the splitting amplitudes, this 
implies that a natural way to extract infrared-finite quantities from non-MHV amplitudes is 
to divide them by the MHV amplitudes with the same number of external legs. Drummond, 
Henn, Korchemsky, and Sokatchev (DHKS) showed to one-loop order that this ratio is not 
only finite but also dual conformal invariant for NMHV amplitudes and conjectured 

that the same holds to all orders for all non-MHV amplitudes [2^. Here we clarify and 
test this conjecture to two-loop order for the six-point amplitude. This test requires the use 
of Dixon and Schabinger's result 25j for the 0(e) terms in the one- loop NMHV six-point 
amplitude. 



As an intriguing consequence of the semiclassical approach of Alday and Maldacena , 
anticipated by the structure of flat space string theory scattering amplitudes at high energy 
and fixed angles 26] , to leading order in the strong coupling expansion all scattering ampli- 
tudes are (in a certain sense) insensitive to the flavor and polarization of external legs. While 
quantum corrections are likely to alter this conclusion, this structure is surprising from the 
standpoint of the intricate analytic structure of the weak-coupling scattering matrix. 

The arguments of Alday and Maldacena led to the identification 16j, [l7j of a surprising 



relation between one-loop MHV amplitudes and the one-loop expectation value of special 



null polygona 



four-, five- 



Wilson loops. This relation was shown to hold at two loops as well for 



and six-particle scattering amplitudes 19 



higher-point two-loop MHV amplitudes are also known 



271 ] . Integral representations of 



; comparison with Wilson-loop 



expectation values [29] is hindered, however, by the complexity of evaluating the required 



higher-point two-loop Feynman integrals. Dual conformal symmetry jl?! 



18 



24 



30 



plays an important role in the relation between MHV scattering amplitudes and Wilson 
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loops. This symmetry is manifest for the integrands of both MHV scattering amplitudes 
and Wilson loops, but it is broken by the dimensional regulator. Dual-conformal invariants 
can be constructed by using the general structure of divergent terms. A particular pattern 
of spontaneous breaking of the gauge group provides an alternative regularization in which 
this symmetry is restored through natural transformations of the regulator [32]. 

It would be interesting to understand whether non-MHV amplitudes also exhibit a similar 
presentation in terms of Wilson loops. A necessary condition is that they exhibit dual 
conformal invariance upon extraction of infrared divergences. It is possible to argue that, to 
all orders in the loop expansion, four-dimensional cuts of any planar scattering amplitude 
in M = 4 SYM, in particular non-MHV amplitudes, have this symmetry. Hints in this 
direction also come from the Grassmannian interpretation of leading singularities; in that 
framework it was shown 33|, |34j that leading singularities are dual conformally invariant. 



Whether this symmetry survives in the complete amplitude, in the presence of the terms 
not constructible from four-dimensional cuts, is an open question. Here we will see that the 
parity-even part of the six-point NMHV amplitudes can be expressed in terms of pseudo- 
conformal integrals, i.e. dimensionally regulated integrals that are invariant under dual 
conformal transformations when continued off-shell. 

While the structure of collinear limits of non-MHV amplitudes is somewhat more intricate 
than those of MHV amplitudes, the former are governed by the same splitting amplitudes 
as the latter. The iteration relation for MHV amplitudes 9j suggests that one can capture 
both the infrared-divergent parts of non-MHV amplitudes, as well as the amplitudes' be- 
havior under collinear limits, via an exponentiation ansatz for all the scalar functions that 
characterize them. This is similar in spirit to the BDS j§| exponentiation ansatz for MHV 
amplitudes. Such an ansatz is not expected to hold to all orders. Departures from it are 
characterized by dual conformal invariant functions which have properties analogous to the 



MHV remainder function 



19 



27]. 



We perform the calculation using the generalized unitarity-based method, employing a 
variety of four-dimensional and D-dimensional cuts to express the amplitude in terms of 
six-point two-loop Feynman integrals. The four-dimensional cuts are evaluated in on-shell 
superspace [35|. This approach automatically takes into account supersymmetry relations 



between different components of cuts and also offers guidance in organizing the calculation. 
We find that the (appropriately defined) parity-even part of the six-point amplitude may 
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be expressed as a sum of pseudo-conformal integrals [30(, in close analogy with the four- 
point amplitude through five loops jsl, [36^ 39 1 and the parity-even part of the five-point 
amplitude through two loops 40M43| . There are some additional integrals in the one- and 
two-loop six-point amplitudes, whose pseudo-conformal nature is less clear. Their integrands 
vanish as D — > 4, yet their integrals can be nonvanishing in this limit. We evaluate the 
integrals using the AMBRE 44J and MB 45j packages and compute the amplitude numerically 
at several kinematic points, related in pairs by dual conformal transformations. The infrared 
singularities of our expression have the structure expected from general considerations j^l 
0] . We have tested numerically the dual- conformal properties of the various finite functions 
that can be constructed from the six-point NMHV amplitude. 

The paper is organized as follows. We review the tree-level and one-loop six-point ampli- 
tudes in section HIl along with their superspace presentation and their conjectured properties. 
Most importantly, we identify a canonical separation of the six-point NMHV amplitude into 
parity-even and parity-odd components. We expect this separation to extend to all orders 
in perturbation theory. In section IHIl we discuss the expected structure of the six-point 
NMHV amplitude to all loop orders, based on our calculation using generalized unitarity. 
We introduce certain finite functions that characterize the amplitude and are expected to 
be invariant under dual conformal transformations. In section IIV} we describe some of the 
details of our calculation. We use a superspace version of the generalized unitarity method. 
We discuss some of the subtle points, and give details on the calculation of two important 
cuts. In section |Vl we present an integral representation of the even part of the two-loop 
six-point NMHV amplitude. For completeness we also list the even part of the two-loop 
six-point MHV amplitude in our notation. We proceed in section IVT1 to analyze our analytic 
and numerical results for the amplitude, and to test the dual conformal-symmetry properties 
of the various functions that have been conjectured to be invariant under dual conformal 



transformations. We give our conclusions and a selection of open problems in section IVIII 
Note added: As the writing of this paper was being completed we received ref. [92] in 
which an alternative presentation of the six-point NMHV amplitude was proposed as a 
consequence of a generalization of the Grassmannian duality for leading singularities to the 
full amplitude. The result also contains a proposal for the parity-odd part of the amplitude. 
Unlike our result, it is expressed in terms of a basis of chiral, tensor integrals written in 
momentum-twistor space. 
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II. REVIEW 



The n-point L-loop planar (leading-color) contributions to scattering amplitudes of an 
SU(N C ) gauge theory with fields in the adjoint representation may be written as 1 



A^ = a L ^[T^...^^^^!)^^);...;^)^^)), (2.1) 

p£S„/Z n 



where we 
refs. 



s 



"ollow the normalization conventions 2 of ref. 



19| (which differ from those used in 



24|). The loop expansion parameter a is, 



(4^rA = (4*0-<^. (2.2) 

^>7l Z 87T 



Here A is the 't Hooft coupling constant and 7 is the Euler constant, 7 = — r'(l). The 
sum runs over all the noncyclic permutations of the external legs, each of which carries 
momentum hi and a polarization vector £j. 

Choosing a specific helicity and flavor configuration for the external legs reduces 
An (k p (i), £ p (i)", ■ ■ ■ ; k p ( n ), £ p (n)) to a color-ordered partial amplitude. Every partial ampli- 
tude can be decomposed into a sum of terms, each of which is a product of a function 
ensuring the correct transformation properties of the amplitude under Lorentz transforma- 
tions (henceforth called "spin factor") and a (pseudo-)scalar function which may be written 
as a sum of L-loop Feynman integrals (the "loop factor"). The spin factor is a rational 
function of the momentum spinors Aj and Aj associated to the external legs; the parity-even 
parts of the loop factor are functions of external Lorentz invariants alone, while the parity- 
odd parts also depend on Levi-Civita contractions of the external momenta. One could of 
course choose to re-express the Levi-Civita contractions in terms of spinor variables. 

MHV amplitudes have two negative-helicity, and any number of positive-helicity, external 
legs. These amplitudes in MSYM have the simplest structure of all amplitudes: they have 
a single spin factor, which is equal to the tree-level scattering amplitude. Computing the 
L-loop MHV amplitude thus amounts to finding the ratio 

^(L),MHV 

M n L) = ,(0),MHV • ( 2 ' 3 ) 



1 We normalize the classical action so that the only coupling constant dependence is an overall factor of 

-2 
9ym- 

2 This definition of the loop expansion parameter extracts the complete dependence on the Euler constant 
from the momentum integrals. 
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CPT implies similar properties for the MHV amplitudes; they also contain a single spin factor 
which is the tree-level amplitude and their scalar and pseudo-scalar factors are obtained from 
corresponding MHV amplitude by a parity transformation. (For alternative presentations 



of the parity-odd terms in M n in terms of spinor variables, see refs. [48|, |49|.) 

All-gluon NMHV amplitudes have three external legs of negative helicity, and any number 
of positive helicity. They are the next-simplest amplitudes after the MHV ones. The five- 



point NMHV amplitudes are MHV; the simplest distinct ones appear for six external legs. 
These have three independent helicity configurations. In contrast to the MHV amplitudes, 
NMHV amplitudes contain several distinct spin factors; their forms depend on the helicity 
configuration of the external legs. As a consequence of relations between spin factors, there 
are many possible presentations of the tree-level amplitudes. We can single out a canonical 
form by constructing the corresponding one-loop amplitude and taking the form that appears 
as the coefficient of the double pole in the dimensional-regularization parameter e. This 

nteed by the general theorems governing the factorization of soft 



relation 



and collinear divergences. We will focus here on the the six-point amplitude. 



A. The Six-Point Gluon Scattering Amplitude at One Loop 

All six-gluon NMHV amplitudes may be obtained by applying cyclic permutations and 

reflections to the three independent helicity configurations (++H ), (+H 1 ) and 

[-\ 1 1 — ). The one- loop amplitudes for these configurations were first obtained in ref. {^J 



through £>(e°) (see also ref. |22j). They can be expressed in terms of three different spin 



factors. The spin factors for the 'split-helicity' configuration (++H ) are, 



r> _ • ^123 (9 A) 

1 (12)(23)(l(2 + 3)4](3(l + 2)6][45][56] ' K ' } 



B =i (4(2 + 3)l] 3 (56)- ^ ( , 

(23) (34) (2(3 + 4)5] [56] [61] s 2M (61) (1(2 + 3)4] (5(3 + 4)2] [34] s 2U ' 1 ' ; 

B _ . (6(1 + 2)3] 3 (45) 3 [12] 3 

3 (61) (12) (2(1 + 6)5] [34] [45] s 3 4 5 (34) (3(1 + 2)6] (5(1 + 6)2] [61] s U5 ' 1 ' ) 

we refer the reader to the original paper |2| for the spin factors of the other independent 
helicity configurations. In all cases, the spin factors are uniquely determined by cuts in 
three-particle invariants. 

The six-point one-loop NMHV amplitude for the (++H ) helicity configuration is 
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jlm.penta 




FIG. 1: The integrals contributing to the six-point one-loop MHV and NMHV amplitudes. An 
arrow marks the leg with momentum fci; the remaining momenta follow clockwise. The one-mass 
box I lm and two-mass easy I 2me integrals contribute to the MHV amplitude and the one-mass box 
I lm and two-mass hard J 2mh integrals contribute to the NMHV amplitude. The one-mass pentagon 
jim,ponta an( j hexagon J hex have numerator factors of (the square of the (— 2e)-dimensional 
components of the loop momentum), and hence are finite. They contribute to both the MHV 
and NMHV amplitudes only at 0(e) and higher (J hex contributes to the even parts while / lm >P enta 
contributes to the odd parts). 



given by, 



where, 



4 1) ' MHV (+++ ) = ~ (bM 1) + B 2 W? + BsWM) + 0(e) 



(2.7) 
(2.8) 



(2.9) 



W i 1] = -\H (l s ^ llm (°) + ~*6i W 2m V)) + 0(e) ; 

and the sum runs over the permutations, 

Si = {(123456), (321654), (456123), (654321)} . 

All permutations in S\ leave the spin factor Bi invariant. The integrals in eq. f)2.8p are 
shown in fig. [lj The factors of ^ in the summand in eq. (12. 8 p are symmetry factors needed 
to compensate for double counting in the summation over Si. The expressions (12.71) and 
( 12~8|) hold only through order 0(e°). At 0(e) eq. (EH]) receives contributions from additional 
integrals while equation (12. 7p receives contributions from additional spin factors. The terms 



25|]. 



of higher order in e have been computed only recently 

The other two scalar functions, and W% \ may be obtained from eq. (12. 8p by re- 
placing the set of permutations S\ by the sets S2 and 1S3, respectively, where 

5 2 = {(234561), (432165), (561234), (165432)} , (2.10) 

5 3 = {(345612), (543216), (612345), (216543)} . (2.11) 



S 



The elements of each of the permutations sets Si, S2 and 1S3 leave invariant the spin factors 
Bi, B 2 and B3, respectively. The union of these three permutations sets, So = S1US2US3, is 
the set of all cyclic permutations and their reflections; the MHV amplitude can be expressed 
as a sum over this larger set. 

The one-loop scattering amplitudes for the other two independent helicity configurations 
have a structure similar to eq. (12.71) : the scalar functions are unchanged while the spin 

factors Bi, B% and B% are replaced 2| by new spin factors Di, D2 and D3 for the (-H 1 ) 

helicity configuration, and by Gi, G2 and G3 for the (H 1 1 — ) configuration 3 : 

4 1),NMHV (++-+— ) = \ + D 2 W 2 W + DM 1] ) , (2-12) 

4 1),NMHV (+-+-+-) = \ (gM 1] + G 2 W? + G 3 W^) • (2-13) 

Infrared consistency then implies that the tree-level amplitudes for the corresponding helicity 
configurations are [2j, 

4° ),NMHV (+++ ) = \ {Bi + B 2 + B 3 ) , (2.14) 

4 0) ' NMHV (++-+— ) = I Pi + D 2 + £3) , (2.15) 



4 0) ' NMHV (+-+-+") = \ (Gi + G 2 + G 3 ) . (2.16) 



2 
1 

2 

The classic expression for these amplitudes was derived in ref. |55j. In later sections we 
will see that the structure present in eqs. (12. 7p . (I2.12p and (I2.13P — in which only the spin 
factors change between various helicity configurations of the external lines — persists at 
higher loops as well. 

B. Superspace and Superamplitudes 

On-shell superspace provides a very convenient way of organizing amplitudes in M = 
4 SYM theory and making manifest supersymmetry relations between them. The bosonic 
part of this superspace is parametrized by the usual bosonic spinor variables Aj, Aj, related to 



3 The notable difference between Bi and the non-split helicity spin factors Di, Gi is that, while the former 
are rational functions of products of adjacent spinors, the latter also contain products of non-adjacent 



spinors. This obscures their transformation properties under the dual conformal symmetry 5J], which 



become manifest only when the amplitudes are combined into a superamplitude 24 
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the external momenta fcj by ki a a = \ia\a- The fermionic part is parametrized by Grassmann 
coordinates rjf, where A = 1, • • • , 4 is an i?-symmetry index. The on-shell fields of the J\f = 4 
theory are assembled into a superfield, 

Hv) =9-+ V A *Pa + ^V A V B (Pab + ^V A V B V C £ABCDiJ D + -^v A V B V C V D tABCDg+ ■ (2.17) 

A superamplitude is a generating function for the scattering amplitudes of component 
fields, which may be identified as the coefficients of the appropriate combinations of r\i 
variables. 

The component amplitudes may be extracted by multiplying the superamplitude with 
the appropriate superfield and integrating over all Grassmann parameters: 

/n 
\\d\ W^hiii) Ai(&l,77l,---,&n,?7n) ■ (2.18) 
i=l 

The superfields $^(7?) have a single nonvanishing term corresponding to the field with 
helicity h^. 

As an example, the n-point NMHV gluon scattering amplitudes appear inside the super- 
amplitude as follows: 

A n (h,T]l,---,kn,Vn) = +VlV2Vs A n(-, ~, ~, +,+,•••, +) (2-19) 

+VlV2V4 A n(-, -,+,-,+,...,+) + •••, 

where r] 4 is the S'f/(4)-invariant expression ^€abcdV A V B V C V D ■ m extracting these compo- 
nent amplitudes, the 77 variables corresponding to the positive-helicity gluons are supplied 
by the superfields (12.171) while those for the negative-helicity ones appear explicitly in the 
superamplitude. Because the half of the supersymmetries manifest in this on-shell super- 
space can be preserved at all stages of scattering amplitude calculations, eq. (12 . 1 9[) holds to 
all orders in perturbation theory. 

The dual superspace in which the superfield is given by, 

®(v) =9+ + Va^ A + ^VAVB<P AB + ^VAVBVce ABCD iJD + ^VAVBVcVD€ ABCD g- , (2.20) 



has also been used, see for example refs. 5|, |6|, |24j, |56| . While the expression for the superam 



plitude is unchanged, component amplitudes are extracted by differentiating with respect 
to selected superspace coordinates, eight for MHV amplitudes, twelve for NMHV ones, etc.: 

A n (h, hi; . . . ; k n , h n ) = $ Aj (Jj^j • • • , K, Vn) ■ (2.21) 
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For pure-gluon amplitudes, the differentiation is solely with respect to the Grassmann coor- 
dinates of the negative-helicity gluons. The structure of superfields is, however, unimportant 
for the computation of superamplitudes. 



In general, n-point tree-level scattering amplitudes can be written as follows 

,(o) _ ^ (4) (Er = i AiAQ^)(Er =1 A^) n ^ vk (9 29] 

(12}(23}-.(nl) ^ 

where V* are polynomials in the Grassmann variables rji of degree 4k. Invariance under 
i?-symmetry implies that are invariant under 577(4) rotations of the Grassmann vari- 
ables r] A . The lowest-order term in the 77 expansion has Grassmann weight 8, while the 
highest-order term has Grassmann weight 4n — 8. CPT conjugation exchanges weight 4k + 8 
with weight 4n — 4k — 8. The k = term in eq. (I2.22p has V® = 1 and contains all the MHV 
amplitudes. The NMHV amplitudes are contained in the k — 1 term. Similarly to equation 
(I2.19p and for the same reason, equation ( 12. 22 ft is expected to hold to all orders in pertur- 
bation theory. Higher-order corrections can alter only the coefficients of the polynomials 
T 3 *, i.e. the component amplitudes. Throughout the paper, four-fold bosonic momentum- 
conserving delta functions will appear, products of delta functions over the four components 
whose indices (a vector index /iora pair (a, a) of spinorial indices) we suppress. A variety of 
four-fold Grassmann delta functions, products of delta functions taken over the SU(4) index 
A, and eight-fold Grassmann delta functions, products of delta functions taken over a pair 
of a spinor index a and an SU(4) index A, will also appear. In these delta functions, we will 
suppress the (boson,c) spinor Mex, but d, S play the (Grassmann) SU(4) index explicitly. 
The tree-level MHV superamplitude was written down long ago by Nair [35|, 

^' = V (12> -(nl> ' (2 ' 23) 



The MHV amplitude has an equally simple form in the conjugate superspace, whose coordi- 
nates are the conjugate spinors A, and the Fourier-conjugate fj of the Grassmann variables 
rj. Fourier-transforming to the same superspace as the MHV amplitude implies [5j that the 



MHV superamplitude is 



^(o),mhv = /^"^^ / d*uf[6W(ri? - A>f). (2.24) 
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Manifestly supersymmetric expressions for non-MHV amplitudes could be obtained 57 1 



through a supersymmetric generalization of the MHV vertex expansion |58M60| . The expres- 
sions obtained this way do not a priori exhibit any special properties. DHKS presented 24 1 
a special form for V$, and showed that it enjoys an extended symmetry, so-called dual 
superconformal symmetry Explicit expressions for all the polynomials were given by 
Drummond and Henn 54J, using a supersymmetric form [6ll . |62| of the Britto, Cachazo, 
Feng, and Witten (BCFW) on-shell recursion relations 

On-shell superspace encodes the relations between amplitudes that are implied by su- 
persymmetry, but does not identify the basic, irreducible components from which all others 
can be obtained. Identifying such basic amplitudes, from which all others can be obtained 
via supersymmetry transformations (along with the required sequence of transformations) 
is in general a difficult problem. Not all corrections to the coefficients in the polynomials 
V* are independent; as these coefficients are nothing but component amplitudes, they are 
related by supersymmetry Ward identities. Elvang, Freedman, and Kiermaier have provided 



a solution [63[ to this class of questions. 

Apart from clarifying the structure of tree-level amplitudes, knowledge of tree-level su- 
peramplitudes allows us to perform manifestly-supercovariant higher-loop calculations using 
generalized unitarity. The one-loop calculation of ref. [24} generalizes the result of ref. [2j for 
the NMHV six-gluon amplitudes to a manifestly supersymmetric expression encompassing 
all possible external states. In refs. Q superamplitudes were used to evaluate the sum 
over all the particles crossing generalized unitarity cuts for n-point MHV amplitudes at any 
loop order. In the section |IV] we will describe in detail the steps needed for evaluating two- 
and higher- loop superamplitudes for any number of external legs and Grassmann weight, 
and elucidate the subtleties that arise in such evaluations. 



C. The six-point NMHV superamplitude 

As our focus in later sections will be on the two- and higher-loop six-point NMHV (su- 
peramplitude, we first review and extend the supersymmetric results of ref. [24) for the 
tree-level and one-loop expressions for this amplitude. 

As is true for the component amplitudes, relations between rational functions of bosonic 
spinor products and Grassmann variables allow the tree-level superamplitude to be expressed 
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in several equivalent forms. We may identify a canonical form, which will also be useful for 
higher-loop calculations, by starting from the e-pole terms in the one-loop superamplitude. 



This superamplitude is given by the supersymmetrization [24| of eqs. (12. 7p . ( I2.12p and (I2.13p . 



4 1),NMHV = ^4 0) ' MHV ((^i3+^^ , 

(2.25) 

where A^' MKY is the tree-level MHV superamplitude, the loop expansion parameter a is 
defined in eq. (12. 2p and the products A^' MHV Rjj +3 j +5 with j — 1, ... ,6 (all indices under- 
stood mod 6) are, 



^(0),MHVd _ 

5 (8) (EW) 



(j (J + 1)> (O'+l) (j+2)> [(j + 3) (7+4)] [(j+4) (j+5)\ 
^ (4) (Vf + s [(J+4) + V f +A [{j+5) (i+3)] + rjf +5 [(j+3) (j+4)]) 



(2.26) 



01 K j+ltj+2 \(j+3)} ((j+2)\K j+3d+i \(j+5)} s jiJ+1J+2 

This product is covariant under dual inversion, with the same weight as the tree-level MHV 
superamplitude. For generic momentum configurations (that is, away from soft or collinear 
configurations), the superfunctions Rjj +3 j +5 are thus invariant under dual superconformal 
t r ansf ormat ions . 

The functions have identical poles in the dimensional regularization parameter e; 
this reflects the universality of infrared divergences. A canonical expression for the tree-level 
six-point NMHV superamplitude is then simply, 

^(o),nmhv _ _^°) ,MHV (jj; 146 -f ^ 251 -f ^ 362 -f R 413 -\- R 524 -|_ R m5 ) . (2.27) 

The R invariants are not all independent; in the presence of the super- momentum conser- 
vation constraint they obey the linear six-term relation, 

*46°^ MHV (-Ri46 — -R251 + -R362 — -R413 + -R524 — -^635) = . (2.28) 
This relation, akin to relations derived from the Grassmannian formulation of tree-level 



amplitudes [65|, leads to two apparently different presentations of the six-point NMHV 



superamplitude: 

^(o),nmhv _ ^°)' MHV ^ 146 _|_ ji m2 _|_ R 524 ) = ^4g°' ) ' MHV (i?25i + -R413 + -^635) • (2.29) 
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A proof of eq. (I2.28|) amounts to showing that the first expression in eq. 



derived from the BCFW recursion relation 



a 



with a supersymmetric shift 



1X29]) can be 
62J while the 



second expression follows from the cyclic symmetry of the superamplitude. At higher loops 
the identity (12.281) is crucial for ensuring the consistency of unitarity cuts. It should also 



play a role in reconstructing scattering amplitudes from their leading singularities 65]. 

The six-point NMHV amplitude is special among NMHV amplitudes as it exhibits a 
discrete invariance related to parity transformations. We will discuss this symmetry and its 
consequences in the following. A similar discussion generalizes to the 2n-point N n_2 MHV 
amplitudes. As mentioned previously, (CPT) conjugation of superamplitudes amounts to 
Fourier-transforming the Grassmann coordinates (reversing the helicities of all component 
fields) and exchanging spinors and conjugate spinors, Aj 4-> Aj. It is easy to check that this 
sequence of transformations maps the products Af^' MHV Rijk into themselves up to a cyclic 
permutation by three units: 

4 0) ' MHV ^i46 4 0) ' MHV ^ 13 , etc. (2.30) 

This is the supersymmetric generalization of an obvious invariance of the six-gluon NMHV 
scattering amplitudes. Invariance of the six-point superamplitude under this transformation 
in turn requires that the functions Wj^ in equation (I2.25P be invariant under conjugation. 

Apart from terms proportional to the sum of R invariants, the 0(e) part of the one- 
loop amplitude also contain terms which are proportional to differences of R invariants. 
They have been computed directly in a one-loop calculation j^J, and their existence may 
also be inferred from the two-loop calculation we will describe in later sections. As such 
differences are odd under conjugation, they must be accompanied by parity-odd (pseudo- 
scalar) functions W i . 



D. Dual Conformal Invariance and the Six- Point Superamplitude 



As mentioned above, DHKS showed 24j that tree- level amplitudes are covariant, with 



weight (—1), under dual superconformal symmetry. This property extends to the rational 
functions A^' MRV Rijk- To what extent does the symmetry extend to the full one-loop 
amplitude? 

The dual conformal and dual superconformal symmetries are only defined in four dimen- 
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sions. One possible extension is the notion of pseudo-conformality: were we to regulate the 
integral functions by off-shell continuation, they would become dual conformal invari- 
ant, as they are sums of box integrals with the appropriate prefactors. Additional evidence 



towards a kind of dual conformal invariance comes from the observation 
singularities are dual conformal invariant. 



33 



34( that leading 



24| that the ratio 4 of the six-point NMHV to 



We can do better than this. DHKS noticed 
MHV superamplitudes, each taken through one-loop order, is invariant under dual conformal 
transformations. That is, the ratio is invariant under transformations that preserve the 
cross-ratios 

S 12<S45 $23 s 56 S34S6I . . 

U 1 = - - , u 2 = - ; — , U 3 = - - . (2.31) 



<Sl23 s 345 S234S123 S345S234 

The ratio of superamplitudes is a natural quantity, as it is infrared finite. 

In gauge theories, the structure of infrared divergences in dimensional regularization is 
independent of the helicity configuration Q, [47]. At one loop, for example, the pole terms 
are proportional to the tree amplitudes. This makes the ratio of any helicity amplitude to 
the MHV amplitude infrared finite. 

The finiteness of such ratios makes it possible to take the four-dimensional limit, and 
to inquire about their properties under dual (super) conformal transformations. Of course, 
finiteness does not guarantee dual conformal invariance. Indeed, the relation between these 
two properties has been investigated in ref. 5| with the conclusion that, in dimensional 
regularization, there exist infrared-finite combinations of pseudo-conformal integrals which 
are not dual conformal invariant. 

Explicit calculations show that such subtleties do not arise here and, through one-loop 
order, the six-point NMHV superamplitude has the factorized form 

1 



24) 



A 



NMHV 



^A 



MHV 
6 



#146 ( 1 + aCS 6 ) + cyclic + 0{af 



(2.32) 



with the functions 6^+3^+5 manifestly expressed in terms of the dual conformal ratios (j2.31j) : 

3 

2 



1 3 

C{% = - In mi In u 2 + - ^2 (In w/c In w fc+ i + Li 2 (l - u k )) - 



k=i 



71 

~6~ 



0(e) 



etc. 



(2.33) 



4 This ratio is sensible because in chiral on-shell superspace any superamplitude is proportional to the super- 
momentum conservation constraint 5^ (X]™=i ^ivf)? which contains the entire Grassmann-dependent fac- 
tor in the MHV amplitude. 
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This function differs from V^ 1 ) as denned in ref. |24J by — tt 2 /6, due to differences in nor- 
malization of amplitudes and finite differences between the Wilson loop expression and the 
one-loop amplitude. It also differs in including 0(e) terms; its e-independent part agrees 
with the function defined in ref. jlj]. 

For completeness we record jl| the integral representation of the one-loop six-point MHV 
amplitude through O(e ): 

^(1),MHV = ^(0),MHV M (1) ; {2M) 



M 6 « 



o-e5iU<s 2 u5 3 

In writing eq. (I2.33j) we used the convention that u i+3 = Ui. The ratio function is thus 
manifestly dual conformal invariant through one loop. It does not have the full dual super- 



conformal invariance, dual supersymmetry being broken by a holomorphic anomaly 
DHKS conjectured 



661. 



241 ] that the main features of eq. f 1 2 . 3 2 j) survive higher-loop corrections: 



that the six-point NMHV superamplitude may be factorized as 

^NMHV = l^MHV ^NMHV + Q ^ . (335) 

and that the functions i?^ MHV have no further e dependence, are well-defined in four di- 
mensions and, to all loop orders, are dual conformal invariant. The conjecture does not 
specify the structure of the 0(e) terms or of the spin factors that enter the functions i?^ MHV 
beyond one-loop level. At one-loop, the 0(e) terms are irrelevant to any 'physical' quantity. 
However, these terms will contribute nontrivially to both the divergent and finite parts of 
the 0(a 2 ) terms in the product on the right-hand side of eq. (12.351) . Our calculation will 
clarify the meaning of these one-loop terms for that part of the amplitude dependent on 
parity-even combinations of R invariants. We will show that they are determined by the 
(9(e) terms in the one- loop NMHV amplitude, which have been calculated recently by Dixon 
and Schabinger j^. 

Before proceeding to describe our calculation, we will discuss in the next section the 
structure of our result as well as the expected properties of the resummed six-point NMHV 
amplitude. 
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I I I I II 

(a) (b) 

Ml Ml 

(c) (d) 
(e) 

FIG. 2: Generalized cuts required to determine the two-loop NMHV amplitude: (a) the 'double- 
pentagon' cut (b) the 'turtle' cut (c) the 'hexabox' cut (d) the 'flying-squirrel' cut (e) the 'rabbit- 
ears' cut. Unlike the MHV calculation, all permutations of the external legs must be considered. 



III. STRUCTURE OF THE SIX-POINT NMHV AMPLITUDE 

In order to obtain the six-point NMHV amplitude to a given loop order, we must de- 
termine all spin factors that appear, and construct the functions of external momenta and 
coupling multiplying each one of them. In the next section we will show explicitly that, 
through two-loop order and through O(e ), the R invariants are the only spin factors that 
appear in the superamplitude. The transformation of the R invariants under conjugation 
( I2.30p then implies that, through two- loop order, the superamplitude can be written as 
follows, 



jNMHV _ 1 yi(0),MHV 
•^6 ~ 2 



(#413 + Ri4s)W 1 (a) + (#524 + R2 5 i)W 2 (a) + (# 635 + R 362 )W 3 (a) 
(#413 - RimW^a) + (# 524 - R 251 )W 2 (a) + (# 635 - R 362 )W 3 (a) + 0(e)! , (3.1) 
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where the Wi(a) are scalar functions and the Wi(a) are pseudoscalar functions. We present 
the calculation of the two-loop six-point NMHV superamplitude, computing explicitly the 
terms depending on parity-even combinations of R invariants. We will find that the four- 

(2) 

dimensional cut-constructible part of the parity-even functions can be expressed as a 
sum of pseudo-conformal integrals. We will also confirm that, unlike their one- loop counter- 

(2) 

parts, the pseudoscalar functions W£ have nonvanishing divergent and finite parts in the e 
expansion. We will not compute these functions explicitly, but the general infrared structure 
of gauge theories divergences requires that they have at most simple (1/e) poles, as both 
the tree and one-loop amplitudes [through O(e )] are free of such terms. In this section, we 
describe the expected general structure of the NMHV amplitude, and the structure of its 
collinear limits. 



A. Beyond Two Loops 

We expect the pseudo-conformality of the coefficient functions to continue to all loop 
orders. To see this, consider a four- dimensional generalized unitarity cut that decomposes 
an L-loop superamplitude into a product of k tree-level superamplitudes 



=T[A 1 A k (3.2) 

cut 

As mentioned earlier and shown in 



24|, each superamplitude has weight (—1) under dual 
inversion. Because a cut propagator simply identifies the Grassmann variables and momenta 
of the legs that are sewn, it has weight (+2) under this transformation. Thus, the product 
above together with the cut propagators has vanishing weight for the cut legs and weight 
(— 1) for the external legs. This implies that these cuts can all be saturated by cuts of 
pseudo-conformal integrals. 

(2) ' — (2) 

Unlike the scalar functions W- , the pseudo-scalar functions W\ are not uniquely de- 
fined. Indeed, the identity ( 12. 28ft implies that it is possible to uniformly add an arbitrary 

' — (2) 

pseudoscalar function to the without affecting the value of the amplitude. In particu- 
lar, we could set any one of these functions to zero. This ambiguity can be partly eliminated 
by requiring that the superamplitude be manifestly invariant under cyclic permutations of 
external legs: 

Wf ] =FW?\, (3.3) 



18 



FIG. 3: A cut of an L-loop six-point amplitude isolating an [L — 2)-loop four-point amplitude with 
no external legs. The cut is proportional to a lone R invariant. 



where P is the operation of permutation to the right by one unit: 

P/[fei, k 2 , fe 3 , fe 4 , fe 5 , fe 6 ] = /[fe 2 , fe 3 , fe 4 , k 5 , fe 6 , fei] . (3.4) 

The corresponding equation for the W\ functions, 

Wf ] = FW?\ , (3.5) 

follows from the symmetry of the superamplitude. 

Requiring cyclic symmetry does not completely fix the ambiguity in the pseudoscalar 

' — (2) 

functions , as parity-odd cyclicly symmetric functions do exist. An example of such a 
function is, 

/ — ^1234/1 + ^2345/2 + £3456 / 3 + ^4561 e 5612 ^6123 h, (3-6) 

where /j are parity-even functions of external momenta kj related by the action of the shift 
operator = P/j_i and e ijmn = e^^k^kjk^k^. 

The generalized-unitarity argument above does not reveal the complete set of spin factors 
that appear at higher loops in the six-point NMHV amplitude. The structure of leading 



singularities suggests [33| that new structures beyond the R invariants of one and two loops 
will be generated at three loops for amplitudes with ten or more external legs, but that no 
new structures will appear beyond that order. It also suggests that no new invariants should 
appear beyond two loops for amplitudes with seven or more, but fewer than ten, external 
legs; and that no new invariants will appear beyond one loop for the six-point amplitude. 

We can, however, argue that the spin factors present at one and two loops will appear 
to all loop orders. As we will see in the next section, all tree-level R invariants appear in 
double two-particle cuts in a channel carrying a three-particle invariant. Such a double cut, 
shown in fig. 12(a), isolates a tree-level four-point amplitude with no external legs attached 
to it. An all-loop generalization of this cut is shown in fig. [3j The well-known property 
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of the four-point amplitude at any loop order, that its spin factor is the same as at tree 
level, implies that this cut will generate exactly the same spin factors as at two loops. This 
argument extends trivially to all higher-loop contributions to the six-point amplitude that 
have double two-particle cuts and isolate four- and five-point amplitudes inside them. It 
can be thought of as a direct superspace generalization of the box substitution rule 
At three loops and beyond, however, it is easy to construct cuts that are outside this class. 
Such cut-based arguments thus cannot rule out spin factors beyond the R invariants seen 
to date. 

Apart from terms containing such new spin factors which may start at three loops, the 
organization of the six-point NMHV amplitude in ( 13. ip holds to all orders in perturbation 
theory. It is therefore interesting to discuss the properties of the parity-even and parity-odd 
functions, 

Wi(a) = 1 + aW t {1) + a 2 W^ 2) + . . . and W^a) = 1 + aW t {1) + a 2 W^ 2) + . . . (3.7) 

in equation (13. ip and of the finite functions 

CVi+3,i+5(a) = 1 + aC^ +3ii+5 + a 2 C? +3i+5 + . . . (3.8) 

that appear in the ratio mhv /t4^ hv . The functions (7^+3^+5(0) will not have definite 
parity. This is due to their relation to linear combinations of functions with differing parity 
properties, to wit (W,(a) ± Wi(a)), as well as to division by the MHV amplitude which does 
not have definite parity properties. For later convenience let us introduce the combinations 

Ci(a) and Cj(a), 

2 (^i+3,i,i+2 + C;,i+3,i+5) = Cj(a) 
2(Ci+3,i,i+2 — Ci,i+3,i+5) = Ci( a ) 

The properties of M 6 (a) together with the universality of infrared divergences implies, that 
through two loops, these functions have definite parity up to corrections that vanish in the 
e ->■ limit? 

5 To guarantee that the functions Ci (a) and Ci (a) have definite parity to all orders in perturbation theory 
it is necessary to divide only by the parity-even part of Me (a). 



Wjja) 
M 6 (a) ' 

Wjja) 
M 6 a) ' 



(3.9) 
(3.10) 
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B. Collinear Limits 



The scalar and pseudo-scalar functions Wi(a) and Wi(a) have specific properties dictated 



by the behavior of the amplitude in collinear limits 

L 



Af\. . . , i\ (z + 1) A ^, - - - ) ^ E E SpJit-A^J (* + l) Ai+1 )4 Z_s) (- ■ ■ , *V • • ) , (3-11) 

A=± s=0 

where k = ki + fcj+i and z is the collinear momentum fraction, fcj ~ 2/5. We can rewrite this 
equation for the all-orders amplitude, 



A e (...,i Xi ,(i + l) 



E SpHt_A(*; i\ (< + l) Al+1 )^ 5 (- • • , fc A , . . . ) • (3.12) 



The properties of C;(a) and C;(a) are more intricate as they involve additional contributions 
from M§(a). 

We will find it easiest to discuss the collinear limits in components. Because Wi and Wi 
do not depend on the precise helicity assignment to the external legs, it suffices to discuss 
the split-helicity configuration. In the three independent collinear limits, the spin factors 
in eqs. (12.4j) . ( 12. 5ft and ( 12. 6 ft behave as follows: 

1||2: B 1>3 -> Split trec (l + ,2+,A;-)4 0) (A; + ,3 + ,4^,5-,6-), 5 2 -> 

5 || 6: 5 lj2 ^ Split t ; ec (5-,6-,A; + )4 0) (l + ,2 + ,3 + ,4-,fc-), fi 3 -> 
3 || 4 : £? 2 , 3 -> Spht t + ree (3+,4- ) A; + )4 0) (l + ,2+,A;-,5-,6-) 

+ Split t _ ree (3 + ,4-,A;~)4 0) (l + ,2 + ,A; + ,5-,6-), B 1 -> (3.13) 

The collinear limits of the parity-odd coefficients 5j, which are contained in the parity- 
odd combinations of R invariants, are similar except that the relative sign between the two 
terms in the 3 || 4 limit is reversed. Combining these limits with the overall behavior of the 
amplitude A3. 12j) . we find, 



l\\2: W 1 + W 3 ^r^(z 
W + W 3 ^ r_(z 

5 || 6 : W x + W 2 ->■ r + (z 
W 1 + W 2 ^r + (z 

3 || 4 : W 2 + W 3 -> r + {z 
W 2 + W 3 ^r + {z 



l + ,2~ 



l+,2" 



1",2- 
5~,6" 
3+,4" 
3 + ,4~ 



E(M 5 MHV (A;,3,4,5,6)) 

0(Mf^(A;,3,4,5,6)) 
E(M 5 MHV (l,2,3,4,fc)) 
0(M 5 MHV (l,2,3,4,fc)) 



(3.14) 



E ( M mhv( 1; 2 , k, 5, 6)) + r_(z; 3+ 4") E(M 5 MHV (1, 
0(Mp(l, 2, fc, 5, 6)) - r_(z; 3+, 4") 0(M 5 MHV (1. 



, 2, fe, 5, 6)) 
,2,fc,5,6)) 



21 



in which, 



/ x / a x Split x (z\i x \ (i + l) Al+1 ) 

r- X {z\ i\ * + 1 Ai+1 = ~ .■ tree/' -\ ( ■ , i -.^ ! 5 3.15 

Split" (z; z Ai , (? + 1) A »+ 1 ) 



M 5 MHV and M 5 MHV are the ratios of the resummed five-point MHV and MHV amplitudes to 
their tree-level counterparts; and E and O denote projection operators onto the parity-even 
and parity-odd components. Functions not explicitly mentioned are unconstrained. 

The collinear properties of the functions C i;i+3 ,i+5 and C i+3 ^^ + 2 can be easily found by 
combining equations (13. 9p . f !3.14j) and the collinear properties of the MHV ratio Mq = 
A^^/Aq ,mhv . In particular, they contain the Levi-Civita tensors necessary to transform, 



for example in the 1 || 2 limit, the MHV five-point amplitude factor into an MHV five-point 
amplitude. 



The iteration relation for the rescaled splitting amplitude (13. 15j) js] suggests another 
natural organization of the functions W{, which is similar in spirit to the BDS ansatz {§]: 

oo 

\nW i = Y, al [fM W i 1) M + C i + R $i + °( e )] ■ ( 3 - 16 ) 
i=i 

The structure of infrared singularities and the collinear behavior require that O(e ) and 
0(e l ) terms in the functions //(e) be the same as for the six-point MHV amplitude. The 
functions R§\, are similar in spirit to the remainder function R$ of the six-point MHV 
amplitude. They are closely related to the functions Ci(a) introduced in equation (13.91) : 

d{a) = exp [ 7jr (a) (w^ - il^)] exp [R^a) - R 6 (a)} + 0(e) , (3.17) 

where 7x(a) is the cusp anomalous dimension and R&-i(a) = Y^i=2 al ^i\-> e ^ c - ^ natural 
consequence of the conjecture that Cj(a) are invariant under dual conformal transformations 
is that the remainder-like functions Rq.^o) are also invariant. We will see that this is indeed 
so. 



The expectation [III, [68] that to leading order in the strong coupling limit, all amplitudes 
with the same number of external legs are identical (or, equivalently, that lim^oo lnCj(a) = 
O(a ) rather than 0(y/a)) predicts a simple relation between the remainder functions Re-i 
and the MHV remainder function R$ to this order. Indeed, using the one-loop relation 

W*> 1) - M (1) = C\ 1] (3.18) 
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and the known value of the strong coupling expansion of the cusp anomaly 
that 

R& — R&A V^A 



69 



72j, it follows 



(3.19) 



with cf ) given in eq. (12 .33 p . Using the numerical results presented in later sections one 
may check that the weak-coupling expansion of the ratio appearing on the left-hand side 
depends on the spin factor labeled by i; it seems therefore that a relation of this type may 
hold only in the strong-coupling limit. 



C. Triple-collinear limits 

Multi-collinear limits provide a richer set of constraints on amplitudes with at least six 
external legs. Unlike the collinear limits discussed in the previous section, they probe the 
detailed structure of the dual-conformal invariant functions unrelated to the infrared struc- 
ture of the amplitude. In the case of the six-point MHV amplitude, they provided a physical 



interpretation of the remainder function 19] . The most detailed limit we can consider with 



six external legs involves three adjacent external momenta becoming collinear, 

k a = Zl P, k b = z 2 P, k c = z 3 P, Zl + z 2 + z 3 = l, 0<Zi<l, P 2 ->• . (3.20) 

Let us understand what such limits imply about the six-point NMHV amplitude and, in 
particular, about the remainder-like functions R^a). 
An L-loop n-point amplitude factorizes as follows IotI ]: 

L 

A^(k 1 ,...,k n ^ 2 ,k n ^,k n ) M- ^^4 L - s )(fc 1; ...,P A ) Split^A^A^A^P). (3.21) 

A=± s=0 

Taking into account parity and reflection symmetries, there are six independent triple- 



collinear splitting amplitudes [19] : 



Split + (A;+A;+fc+;P), (3.22) 
Sp]it_ Xp {k^k^ b k^;P), A a + \ b + A c - A P = 2 , (3.23) 

Sp\it_ Xp (k^k^k^; P), A a + X b + A c - X P = . (3.24) 

The first one (I3.22j) vanishes in any supersymmetric theory. The three triple-collinear split- 
ting amplitudes of the second type (13.231) . an example of which is A a = A^ = A c = Ap = 1, 

23 



appear in limits of MHV amplitudes. The TV = 4 supersymmetry Ward identities for MHV 
amplitudes imply that their rescaled forms 6 are all equal, 

Split«(^ c +;P) _ Split? P) (0 s 



Split*? {ktKki; P=f) Splitf P) 5 Sa6c ' Sa6c ' 

These splitting amplitudes are relevant only for NMHV amplitudes with at least seven 
external legs. They do not arise in the factorization of six-point amplitudes, because the 
four-point amplitude entering the factorization (13.2ip vanishes identically. 

The two splitting amplitudes of the third kind f !3.24j) arise only in limits of NMHV 
amplitudes and do not have a simple factorized form similar to f l3.25j) 7 . They are however 
the only splitting amplitudes that can appear in the triple-collinear limit of the six-point 
NMHV amplitude. 

As is true for the tree-level NMHV amplitudes, the splitting amplitudes f)3.24p have several 
different presentations related by potentially nontrivial spinor identities. A canonical one, 
that is useful for our purpose, is obtained from the triple-collinear limit of the six-point 
tree- level amplitude in equations (I2.14p - fl2.16p . 

As the functions Wi are independent of the helicity assignment of the external legs, we 
again discuss only the split helicity configuration. Up to conjugation and relabeling the only 
non-trivial limit is 2 || 3 || 4. With the momentum fractions fc 2 = ZiP, k s = z 2 P, &4 = z^P, 
the spin factors Bi become: 

, Pi (1-^s) 2 
h = -7*r, — i-» 



4 0) (l+P+5-6-) v^I^(23)(yiT[24] + v^[34]) 



B 2 (xAT(24) + v^<34» 



3 



4 0) (l+P+5-6-) s 234 (23) (34) (V^(23) + ^(24)) 



_j_ H fQ 26) 

s 2 34[34]( v ^r[24] + y^[34])( v ^[23] + v ^[24]) ^ ; 



3/2 



4 0) (l+P+5-6-) V^i 1 ~ ^i)[34]( v ^(23) + ^(24)) 

(^ 3 ) 3/2 



+ 



i 2 -(34)( v ^[23] + v ^[24]) 



We omit a trivial dimensional dependence on s a bc from the argument list of r s . 
7 The spin-averaged absolute values squared of tree-level triple-collinear splitting amplitudes were computed 
in ref. [z3|; without spin- aver aging, in refs. 7_4[. The tree- level triple (and higher) collinear splitting 



a 



mplitudes themselves were computed in ref. [75j using the MHV rules [76|. The one-loop correction to 



the q —> qQQ triple-collinear splitting amplitude in QCD was computed in ref. [77 1. 
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The tree-level splitting amplitude is simply 

Sp\it m (ktktk^;P) = ^(b 1 + b 2 + b 3 ) . (3.27) 

Thus, while these splitting amplitudes do not have a simple factorized form similar to that 
for splitting amplitudes of the second type (|3.25p . we see that the structure of the six-point 
amplitude (13. ip implies that to this order each component 6, is dressed at higher loops by 
scalar functions of momenta, 

Split_(A;+A;+A;4 ; P) = ~ (hw^a) + b 2 w 2 (a) + b 3 w 3 (a)) . (3.28) 

The parity-odd spin factors Bi also have nontrivial triple-collinear limits. Their coefficients 
Wi, though, must contain Levi-Civita tensors and thus naively vanish in this limit. The 
triple-collinear limits of additional spin factors that may appear beyond two loop order 
must be considered separately. 



19j, none of the conformal cross-ratios (I2.3ip 



As was true for the limit discussed in ref. 
vanish as 2 || 3 || 4; they become 

z l z 3 - S 23 1 _ S34 1 , . 

(1 - Zi)(l - Z 3 ) S 234 l~Z 3 S 2U l-Z! 

Thanks to their expected dual conformal invariance (which we will confirm in later sections), 
the remainder-like functions R^i(a) retain their complete kinematic content, and may be 
read off the two-loop triple-collinear splitting amplitude f )3.24p by subtracting the triple- 
collinear limit of the two-loop iteration of the one-loop functions . 

IV. CONSTRUCTING THE EVEN PART OF THE TWO-LOOP AMPLITUDE 

We will construct the even part of the two-loop six-point NMHV amplitude using a super- 



21 



22 



78j. On general grounds, 



space form [5j, |6| of the generalized unitarity method li- 
the result will be expressed as a sum of planar two-loop Feynman integrals with coeffi- 
cients that are rational functions of the spinor variables. At this order, one-loop calculations 

suggest that it is possible to exclude integrals with triangle or bubble subintegrals. 

(2) ' — (2) 

Similarly to the two-loop MHV amplitude, we will find that neither W{ nor W{ can be 
completely determined by four-dimensional cuts. Rather, they receive both divergent and 
finite nontrivial contributions from integrals whose integrand is proportional to the (— 2e) 
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components of the loop momenta. It is quite nontrivial that these latter contributions can 
be organized in terms of the same R invariants as the four- dimensional cut-constructible 
terms. 

The generalized cuts that determine the amplitude are then the ones shown in fig. 121 
which are the same ones that determine the MHV amplitude [19|. Unlike the calculation of 
the MHV amplitude, however, here it is necessary to evaluate cuts with all external helicity 
configurations, as each yields information about different spin factors. 

In any supersymmetric theory the improved power-counting ensures that at one-loop 
order and through O(e ) all terms can be detected in four- dimensional cuts. Beyond one 
loop this is no longer true generically; for example, the six-point MHV amplitude at two 
loops receives nontrivial contributions from integrals whose integrand vanishes identically 
when evaluated in four dimensions. Four-point amplitudes in the M = 4 SYM theory are an 
exception: through five loops they appear to be determined solely by four- dimensional cuts. 
We therefore decompose the functions W> in eqs. ( 13. lj) and ( 13. 71) into a four-dimensional 
cut-constructible part and a part that requires D-dimensional calculations, 

For the former, powerful helicity and supersymmetry methods can be employed. The latter 
part of the amplitude is determined by comparing the result of D- dimensional and four- 
dimensional calculations and is expressed in terms of "/i-integrals" — nontrivial integrals 
whose integrand vanishes identically in four dimensions. 

While all cuts may be evaluated easily, separating their contributions to each one of the 
functions W^' ^ is not always straightforward. As mentioned previously, (multiple) cuts 
in channels carrying three-particle invariants capture a single even and odd spin factor at a 
time and thus determine terms in a single W^' ^ and w[ 2 ' > ' D ~ 4i , with the index i deter- 
mined by the helicity configuration of external legs. This is the case for cuts (a) and (b) 
in fig. [2j In contrast, (multiple) cuts in channels carrying only two-particle invariants con- 
tribute simultaneously to several spin structures and thus to several W^' D ~ A and W^' D ~ 4 
functions. This feature is already present in the cut construction of the one-loop amplitude; 
in that case however, cuts in channels carrying three-particle invariants suffice to completely 
determine the amplitude Q|. The similarity between the expression for cut (c) of fig. [2] and 
a cut of the one-loop amplitude makes it possible to disentangle it. The cuts of fig. Wi.d) and 
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(e) however seem intractable in a component approach. 

The component approach also fails to incorporate in a transparent way the constraints 
imposed by supersymmetry. On-shell superspace provides the additional structure necessary 
for identifying the contributions of the remaining cuts to each of the W- 2S) ' D ~ 4 We shall 
therefore formulate the entire calculation of the four- dimensional cut-constructible part of 
the amplitude in on-shell superspace. After a brief overview of the structure of supercuts and 
of the techniques necessary to disentangle them, we will discuss cut (a), and then proceed 
to a more detailed analysis of the challenging cuts (c), (d) and {e\ 



shall use a superspace generalization of the maximal cut method 
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or the latter cuts we 



79|. 



A. Unitarity in Superspace: General Features and Techniques 

Generalized cuts may be classified following the number of cut conditions they impose. 
The same is true for generalized supercuts. At L loops in four dimensions it is possible 
to impose at most 4L cut conditions; based on one- and two-loop information, it it likely 
that their solutions generically form a discrete set. This type of cut has been considered 
in the maximal-unitarity approach as well as in the leading-singularity approach. Maximal 
cuts, i.e. cuts with the maximal number of cut propagators, are typically insufficient to 
completely determine an amplitude. For example, at two loops one frequently encounters 
double box integrals, which cannot be detected by cutting eight propagators. Near-maximal 
cuts, obtained by successively relaxing cut condition in maximal cuts, provide an algorithmic 
way of identifying these contributions. Near-maximal cuts exhibit additional propagator- 
like singularities which are exploited in the leading singularity approach to reduce the one- 
parameter family of solutions to the cut conditions to a discrete set. 

The two-loop six-point NMHV amplitude can in principle be determined entirely from 
the iterated two-particle cuts shown in fig. The Feynman integrals that contribute only 
to cuts (c), (d) and (e) are also detected by certain near- maximal cuts. We have used them 
instead to check that the resulting amplitude correctly reproduces cuts (c), (d) and (e), 
supplemented by an additional cut condition isolating terms in one of the tree amplitudes. 

General supercuts are constructed 0, [f3, Q by multiplying together superamplitudes, 
identifying the r\ parameters of the lines that are sewn together and integrating over the 
common values of the internal rj variables. The structure and properties of general supercuts 
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have been analyzed in detail in ref. [6| where it was shown that, upon use of a supersymmetric 



generalization of the MHV vertex rules 76), their building blocks are generalized supercuts 
constructed only out of MHV and MHV tree-level amplitudes. 

When evaluating a supercut one encounters the situation that on one side of the cut a 
momentum is outgoing and on the other side it is incoming. In order to write the tree-level 
amplitude and in particular the argument of their delta functions, it is necessary to define 
the spinors \—p) and \—p] corresponding to the incoming momentum (— p). We use the 



analytic continuation rule 56) that the change in sign of the momentum is realized by a 



change of sign of the holomorphic spinor 

p i — y —p -H- A p h> — A p , X p h-> +X P ; 

\-p)*->-\p), \-p]^\p\- (4.2) 

Let us discuss in detail the building blocks we require, supercuts constructed only out of 



MHV and MHV tree-level amplitudes. Their evaluation requires the evaluation of integrals 
of products of delta functions with arguments linear in Grassmann parameters, see eqs. ( 12. 23ft 
and (12.24jl . For a p-particle cut of an N fc MHV amplitude this product contains (8 + 4(/c + »V) 
delta function factors of which (8 + 4k) remain upon integration. As discussed in refs. j^, |fj 
56) . the integration over the internal r\ parameters realizes the sum over the states crossing 
the (generalized) supercut. For any supercut, eight of these delta functions can always be 
singled out: they enforce the super-momentum conservation of the amplitude, 

*< 8 >(X>fl (4.3) 

where £ denotes the set of external lines. These delta functions may be thought of as 
the supersymmetric generalization of the usual momentum conservation constraint. They 
may be extracted without carrying out any Grassmann integrations, by taking suitable 
linear combinations of the arguments of all delta functions. If the supercut contains at 
least one MHV superamplitude factor, the Jacobian of this transformation is unity. Their 
extraction also makes manifest the invariance of the amplitude under half of the maximal 
supersymmetry. Invariance under the other half of the supersymmetry, generated by 

n f) 

is not manifest, but can in principle be checked at the level of the Grassmann integrand. 
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The delta functions (14. 3 j) represent the complete Grassmann parameter dependence of 
a supercut of an MHV amplitude. The Grassmann integrals simply yield the determinant 
of the system of linear equations which are the arguments of the other 4p delta functions, 
where p is the number of cut lines 6]. 

For cuts of an N fc MHV amplitude there is a certain amount of freedom in evaluating the 
internal Grassmann integrals. In general, however, the resulting 4k delta functions have 
many undesirable features. The essential ones are that (1) their arguments may depend 
on loop momenta (if the cut conditions do not completely freeze the momentum integrals) 
and (2) they may not make the symmetries of the amplitude manifest. We wish to express 
these Grassmann delta functions in terms of structures that appeared at lower-loop order; in 
the case of the six-point NMHV amplitude; these are the dual superconformal R invariants. 
This is a non-trivial operation, and we have but a limited set of tools available. 

Given a set of 4k delta functions 

k 

]p (4) (e^,A)), (4.5) 

i=l 

it may be possible to construct linear combinations of their arguments M i j(X)ej(rj, A) which 
factorize into products of the desired combinations of spinors and Grassmann variables upon 
use of momentum and super-momentum conservation, cut conditions, and the fact that a 
Grassmann delta function equals its argument. For k = 1, which is the case of interest to 
us, no linear combinations can be constructed. 

A possible strategy for eliminating the dependence of the delta functions on loop momenta 
is to make use of the fact that a Grassmann delta function equals its argument. This 
observation replaces a cut carrying a Grassmann delta function with a sum of cuts of tensor 
integrals with Grassmann-valued coefficients. Albeit nontrivial due to their high rank, the 
tensor integrals may then be reduced following the standard strategy of integral reduction. 
While indeed successful in eliminating the loop momentum dependence from the Grassmann 
delta functions, this strategy is likely to lead to rather unwieldy expressions. We will not 
pursue this direction. 

An alternate approach to reorganizing Grassmann delta functions is to use the Lagrange 
interpolation formula, which is most efficient when applied to next-to-maximal cuts, which 
impose (4L — 1) on-shell conditions. Let y be the variable that parametrizes the solution 
to these cut conditions. The product of the 4k Grassmann delta functions is then just a 
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polynomial Pd(y) of degree d = 4k with Grassmann- valued coefficients. Any such polynomial 
may be written as 

d+l d+l 

where the values are arbitrary. This equation simply encodes the fact that a polynomial 
of degree d is determined by its values at d + 1 points. 

Choosing the points yi can be regarded as freezing the momentum component unfixed 
by the cut condition; from this perspective it is akin to the leading-singularity method 
which uses additional cut-like conditions for the same purpose. The Lagrange interpolation 
formula (14. 6[) provides a different strategy, as the points yi need not be chosen following 
the leading-singularity prescription. If the two approaches are to agree, the residue of the 
leading singularity must be proportional to the Grassmann delta functions appearing in dual 



superconformal invariants. Evidence that this is indeed true has been presented in ref. (33]. 
In general, however, in order to use the interpolation formula (14. 6p . there must exist more 
Hi such that Pd{Hi) is (proportional to) a dual superconformal invariant than are given by 
the leading-singularity approach. 

In the next subsection we will use this strategy to analyze certain seven-particle cuts 
of the six-point two-loop NMHV superamplitude. As we will see, with judiciously chosen 
points i/i it is possible to have P^iyi) be proportional to the delta functions appearing in the 
R dual superconformal invariants. 

Because of the arbitrariness in the choice of the i/i, the decomposition of Pd{y) in a linear 
combination of "good" Grassmann delta functions, such as the delta functions appearing 
in the dual superconformal invariants, is not unique. This signals the existence of linear 
relations between the dual superconformal invariants. For six-point amplitudes, an identity 
arising this way is eq. (12. 28 h . which was already required for the consistency of the various 
possible presentations of the tree-level amplitude. It is conceivable that at higher points 
and/or higher loops, new relations arise, beyond those that can be obtained from tree-level 
considerations. 
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FIG. 4: Two-loop topologies entering the 2-loop 6-point amplitudes. The arrow on the external 
line indicates leg number 1. 

B. Supercut Example: the Double-Pentagon Cut 

Let us illustrate the general strategy outlined in the previous section, by examining in 
some detail two cuts that are essential for the construction of the six-point NMHV super- 
amplitude. We begin with the 'double-pentagon' cut, shown in fig. E(a), which isolates the 
double-pentagon integrals 1^ and (shown in fig. H]) from a wide class of other integrals 
(hence its name). 

This cut provides two distinct contributions to the coefficients of the NMHV amplitude, 
depending on which of the two five-point tree-level factors is an MHV or an MHV super- 
amplitude. Each of the contributions is closed under supersymmetry transformations, so we 
will call them supersectors. (They were called "holomorphicity configurations" in ref. jf]].) 
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6 1 

(i) (") 
FIG. 5: The two contributions to the 'double-pentagon' supercut[2ja). The circled + and — denote 



MHV and MHV superamplitudes, respectively. The middle amplitude may be chosen to be either 



of MHV or of MHV type. Here we choose to present it as an MHV superamplitude. 



These two supersectors are shown in fig. [5} their values are, 



and 



where 
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(4.7) 



(4.8) 



<zm = E A *v and <=E A ^- ( 4 - 9 ) 

i=l i=A 

By taking appropriate linear combinations of the arguments of the delta functions in these 
equations it is easy to extract the overall super-momentum conservation delta function, 
^ (?123 + ?456)- Carrying out the Grassmann integrals we then find 

Sl23 (l|2 + 3|4](3|4 + 5|6](/ 1 / 2 ) 4 



dp 



5(h) 



A^ MHV Rue 



((3/ 1 )(/ 1 / 2 )(/ 2 i))((/ 2 / 1 )(/ 1 / 4 )(y 3 )(y 2 )) (fihWMM) 

gi 23 (4|2 + 3|l](6|4 + 5|3] (Z 4 Z 3 ) 4 

((6/ 3 ) (hh) (Z 4 4» ((hh) (hU) (kh) (hh)) ([3/i] [hh] [Z 2 1]) 



, (4.10) 
• (4-11) 
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We can reorganize the contributions into even and odd components, 

C d P = G|P + G|P = ^(O)-MHV ^dp ( ^ i3 + R ^ + ^(0),MHV ^dp (jR4is _ Ruq) (412) 



The two functions of vanishing weight in equations (I4.10p and (14. lip may be identified 
as the contribution of gluon intermediate states in a component approach. They can be de- 
composed by standard means, by reconstructing propagators and organizing the numerator 
into a single trace? For example, 



c d P = si23 ( (l|2 + 3|4](3|4 + 5|6](/ 1 / 2 ) 4 

+ " {(hh){hk){Uh){hh)) V«3/i>(W2>(/2l» ([6Z 3 ][y 4 ]M]) 

(4|2 + 3|l](6|4 + 5|3](y 3 ) 4 



1 " 



4 



5l2 3 (5l23 s 345 ~~ S12S45) S 61 S 12 3 



L(/i + h) 2 (h + / 3 ) 2 (/ 3 + h) 2 (h + h)\l 2 + l 3 f(l 3 + h] 

S 34-S 1 23 5i23(Si23S234 — S23 s 5g) 
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(h + h) 2 (l 2 + l z f{h + (l 2 + k x )\k + h) 2 (h + hy 

s 12 s 23 s 123(^6 — ^2) Si2S23Sl23(^4 — h 
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(h + h) 2 (h + k 3 f(h + / 3 ) 2 (/ 3 + h) 2 (h + hnh + k 3 f{l 2 + hf{h + h) 2 

S45S5eSl23(h - U) 2 S 4 5S 5 6Sl23(^l ~ h) 2 



(h + h) 2 (l 2 + h) 2 (h + h) 2 (l 3 + (i 2 + h) 2 (l 2 + h) 2 (U + h) 2 (h + hf 

1 

f- 



(h + h) 2 (h + k 3 f(l 2 + hf{h + h) 2 (h + h) 2 

xsi2 3 (0i 2 S45 _ s i23S345)(^i - hf (h - hf + s 3 iSi 23 (k 1 - l 3 ) 2 (k 6 - hf 
+SQiSi 23 {h - h) 2 (h - hf + (s 23 s 56 - S123S234X&3 - U) 2 (h - h) 2 



+ (k + hnh + k 3 y(u + h) 2 (i 3 + k,y (413) 



X 



(2S12S23S45S56 - Si23(s6lS 3 4Sl23 + S12S45S234 + S23S56S345 ~ S123S234S345 



From this expression we can easily read off the coefficients of all integrals in fig. H] that have 
a double cut in the S123 channel. Some integrals appear multiple times, corresponding to 
different cyclic permutations of external legs that have such a cut. The numerator factors 
in the expression above are precisely those required to render the integrals invariant under 
dual inversion. As we did not need to specify the helicity labels of the external legs, all cuts 
with this topology can be obtained by simple cyclic relabeling. 



This last step is important to avoid the appearance of parity-even terms which are a product of two 
parity-odd factors. 
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(i) (ii) (iii) 

FIG. 6: Next-to-maximal cuts that detect the integrals not easily isolated by the iterated two- 
particle cuts: from left to right, next-to-maximal cuts for the 'hexabox,' 'flying-squirrel,' and 
'rabbit-ears' cuts of fig. [2(c), (d), and (e), respectively. 
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(i) (ii) (iii) (iv) 

FIG. 7: The four possible assignments of internal helicities for the next-to- maximal 'flying-squirrel' 



cut of fig. 6(h) The '0' vertices denote three-point MHV amplitudes while the '©' vertices denote 
three-point MHV amplitudes. 

The 'turtle' cut shown in fig. W(b) can be computed in a similar way, and also contributes 
a lone R invariant. These two cuts determine the coefficients of all integral topologies in 
fig.! except J (7) , J (14) , and J (15) . Other cuts are necessary to determine these contributions. 
An efficient strategy, which makes use of the results obtained from the double-pentagon cut 
of fig. 121(a) and the turtle cut of fig. W[b), is to analyze the relevant next-to- maximal cuts 
and find the remaining integrals one at a time. 



C. Supercut Example: A Contribution to the Flying-Squirrel Cut 

Following this strategy, we present one contribution to the next-to-maximal cut, shown in 



fig. 6 (ii) , that imposes additional cut constraints beyond the 'flying- squirrel' cut of fig. \2^d). 
It serves to isolate one of our target integrals, I^ 7 \ and allows us to determine its coefficient. 
As explained above, we impose the additional cut conditions because of difficulties in or- 
ganizing the results for cuts like the 'flying-squirrel' cut in terms of dual superconformal R 
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invariants. A superspace calculation, combined with the reduced set of Feynman integrals 
isolated by the next-to-maximal cut conditions, reduces the ambiguity in such a reorganiza- 
tion by enforcing supersymmetry relations between the reduced number of contributions to 
these cuts. 



The next-to-maximal 'flying-squirrel' cut of fig. 6(ii) has four supersectors, corresponding 



to three-point amplitudes at the corners being of MHV or MHV type. The supersectors 



are shown in fig. [7J We discuss in detail the configuration in fig. 7(ii) , quote the result for 



the configuration in fig. 7(i) , and explain how to construct the other two components by 
relabeling. 

is 



a 



The product of the tree superamplitudes entering the supersector shown in fig. 7(if 



d r] h d rj h d rj h d r] k d VqA ri q2 d r] m 

[Ui]mh][hM 



(q 1 2)(2l 1 )(l 1 q 1 ) (1,3) (3Z 4 ) (hq 2 ) {q-ih) 

6 {4) (Mv£ + Wtot + Wit) t {8) (Kv£ + W + Kvt 



[l 4 A][Aq 3 ][q 3 l 4 



x 



{q 3 5)(5l 3 }{l 3 q 3 } 



(4.14) 



<6Z 2 ></ 2 g 2 ><^3><«36> 

The expected overall super-momentum conservation constraint may be extracted by adding 
the arguments of all the eightfold delta functions 5^ to the last such function, and then 
using momentum conservation to eliminate A^r/^ and A/ 4 r/^. These transformations have 
unit Jacobian. 

The remaining Grassmann integrals can be computed easily; we obtain: 



Af= [Igi] 4 5 (4) ((gi2)(g 2 / 1 )(g 3 /3)[4g3]^ + (<Mi>(fc3><<feJ 8 )[4<&]rf 
+(qih)(q2k)(q3h)[q3k}vi + (Qik) {qzh) (^5)[/ 4 4]r/ 5 A ) 



(4.15) 



The contribution from the supersector in fig. 7(ii) is obtained by dividing M by the explicit 
denominators in equation (14.141) . This expression can be simplified in several different ways; 
we proceed by solving the cut conditions. The internal spinors (except for those associated 
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(4.16) 



to g 2 ) may be expressed conveniently in terms of two variables y and z: 

K = 2/Ai - A 2 , X h = A 2 , A 9l = yX 1 , 

Az 2 = > A/ 2 = Ai + yX 2 , A gi = A2 , 

A« 3 = —zX^ — A5 , A/ 3 = A5 , X qs = zX^ , 

A; 4 = A4 , A/ 4 = — A4 + zX§ , X q[i = A5 . 

The momentum g 2 can be determined through momentum conservation; the condition that 
it be on shell relates the two parameters y and z. These relations imply that all spinor 
products in eq. (14.151) that do not contain the holomorphic spinor |g 2 ) are monomials in y, 
z and spinor products of external momenta. The remaining holomorphic spinor products, 
which do contain \q 2 ), can be converted into functions of y and external spinor products by 
multiplying and dividing by [g 2 5] 4 and using the identities, 

= -(2(3 + 4)5] +i/(l(3 + 4)5], 
(WM = (16) [65], (4.17) 

(WfeS] = (16)([16]+t/[26]), 
(kq 2 )[q 2 5] = (4(2 + 3)5] +|/(14} [25]. 

The numerator factor M in eq. (14.151) is then, 

N = ^y f(4) (^((2( x + 6 ) 5 ] - y^ 2 + 6 ) 5 ]) + % A «3(i + 6 ) 5 ] - y^)m 

+ r/ 4 A ((4(l + 6)5] - 2/(14) [25]) + % A (s 234 - 2/(1(3 + 4) 2])) . (4.18) 

This expression is invariant, though not manifestly, under the action of the supersymmetry 
generators q = ^ Xd v . 

Overall super-momentum conservation provides the means to further simplify M . By 
subtracting Y^iVii^O- + 6)5] = 0, and adding Y^iWi = to the argument of the 

delta function we find 

(4.19) 

This superspace expression has the two unwanted features already mentioned in sec- 
tion IIV At on the one hand, the argument of the delta function depends on the internal 
momenta through the variable y; on the other, it is not manifestly a function only of the 
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same superspace structures as the tree-level amplitude (13.11) . We wish to reorganize it in 
terms of R invariants, and at the same time remove the dependence on internal momenta 
from the arguments of the delta function by using the Lagrange interpolation formula (14. 6 p 
on the degree four polynomial, 

Pdv) = 5 (4 V[56] + v £[61] + ^[15] + y(t^[56] + r]£[62) + ^[25])) 

5 5 

= EII^%)- (4-20) 

For this to be possible, as explained earlier it is necessary that there exist at least five values 
Ui such that Piijji) is proportional to an R invariant. It turns out that there are at least six 
such values: 

P <(Str||) = ((f^]) 4 ^^[«] + ,f[ 5 3] + ^[34])^3 5 , (4.21a) 

Pi (~w]) = (m ) 4 ' 5(4,(%/l1121 + " i4|2e| + v ° m K R ™ • (4 - 21b) 
SfSS) - ( <fik) V, <^i + + ^ 2 » « *» , <**«> 

A (||) = (H) 4 <5 <4, ta 4 [56] +,, 5 4 [64] + % "[45]) <x fl 146 , (4.21d) 

P * («S4l) ^ ((lfrll)^ (4, (* A i 34 l + nil 42] + rf M « « 5 , , (4.21.) 

P 4 (0) = «**'Cfe A [61] + ^[15] + >)f [56]) cx fl 251 . (4.21f) 

For some of these cases, we have used overall super-momentum conservation constraint as 
well as nontrivial spinor identities to transform the argument of 5^ 4 \ As we will see shortly, 
only the first four values of y correspond to leading singularities. 

We can use the Lagrange interpolation formula for any five of the six special values 

/ (2(3 + 4)5] [16] (4(5 + 6)1] (23) s 234 \ , , 
1^1/^3,^1/5,^- j (1(3 + 4)5] , [26] , ( 4 ( 5 + 6 )2]' (13)' (1(3 + 4)2]' J ' [ } 

Clearly, the decomposition obtained this way is not unique as there are six different possi- 
bilities. Let us denote them by d, where i is the index of the missing root. In general we 
can construct a five-parameter decomposition 

6 

P 4 (y) = ^2a i C i (y), (4-23) 
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with Yl a i = 1- The remaining parameters a, may be constrained by requiring that the 
superamplitude have additional manifest symmetries; for example, one may require that the 
parity of the superamplitude be manifest. We impose such a requirement in the following. 

All the presentations of the cut obtained for different possible choices of five values y,i are 
physically equivalent. However, they contain different R invariants; the existence of more 
than five values is equivalent to the existence of nontrivial relations between R invari- 
ants. These relations hold only in the presence of the overall super-momentum conservation 
constraint. 



We are now in position to assemble the result for the supersector of fig. |7(ii) 



Further use of the identities ( I4.17P implies that it is given by 

/(OJ.MHV D / \ 

-Aq s 12 s i5 P 4 (y) 

(34) [56] (2/(13) - (23))(y<l(3 + 4)5] - (2(3 + 4)5])((4(5 + 6)1] + </(4(5 + 6)2])([16] + </[26]) ' 

(4.24) 

Inspecting the denominator of this expression, we see that the first four points yi in eq. (I4.22j) 
correspond to poles. They are in fact positions of leading singularities, as all of them arise 
from the [(?25]~ 4 factor in eq. (I4.19p which is the Jacobian arising from solving the cut 
conditions. 

Choosing the five points yi to be 2/2, 2/3, 2/4, 2/s} ? we obtain, 

n , rf(o),MHv ( 2/4(2/- vs) „ . 2/2(2/ — 2/5) D 



y 3 (y-ys) ^ _ fttozJ*)^ + ^ , (4 . 25) 



2/5(2/ -2/s) 2/5(2/ -2/1; 

Each of the denominators appearing in this expression may be identified with a propagator 
evaluated on the kinematic configuration ( I4.16p . 

Thus, the contribution of this supersector depends only on R invariants. We can decom- 
pose it in even and odd invariants (-R^+3^+5 ± -Ri+3,i,i+2) , following the form ( 13. lj) of the 
superamplitude. To identify the part of ^7(^)1 that receives contributions from the missing 



integral 1^ we need to subtract from it the contribution of all the other integrals in fig. HJ 
determined from the cuts of fig. |2^a) and (b). These cuts can contribute only terms propor- 
tional to the invariants -R146, -R362, -R413 or i?635- Thus, we can conclude immediately that 
-R524 arises solely as a coefficient of I^ 7 \ whose coefficient must therefore be, 

1 ^( 0)i MHV Sl2S45S234(jR251 + r 52a) . (4.26) 
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Indeed, it is intuitively clear that because of its topology, can appear in the coeffi- 
cient (14.251) only in terms that have no additional propagators. 

Carefully repeating this analysis for the other even invariants implies that the complete 
contribution of this supersector to the even part of /( 7 )'s coefficient is, 

1 „(0),MHV/ D D \ / \ 

— 7*^6 1-0-146 + -n-413j s 123l s 234'S345 ~~ S61S34J 

1 ,(0),MHV/ D d \ / \ 

~4 6 \ n 362 ~r -K635j s 345l s 123 s 234 — S 23 S 5e ) 

+ l^(o),MH V(i?25i + ^ 524 ) Sl2S45S234 . (4.27) 



This conclusion must be checked against the other configurations in fig. [7J Fig. 7(iv 



is 



the parity conjugate of fig. 7(ii) and should therefore yield the same result for the even R 



invariants (and its negative for the odd-parity ones) 



The configurations in figs. 7(i) and 7(iii) and are parity conjugates of each other. Evalu- 



ating them following the same steps yields, 

n ,_ „(o),mhv s 12 s 45 5 (4) ([61k A + l^K + [56]^) i(0),MHVn ( A 0Q \ 

(23)(34)[56][61](4(2 + 3)l](2(l + 6)5] " A R ™ s ^ s ™ ,(4.28) 



n ,_ ,(o),mhv s 12 s^\[U}^ + [42] % A + [23]Q _ 

i(0),MHV p ( A on \ 

- A (16)(65)[23][34](5(l + 6)2](l(2 + 3)4] " A R^s 12 s i5 s 2U .(4.29) 



Unlike the configuration in fig. 7(ii) , the loop-momentum dependence here cancels completely 
after integration over the internal Grassmann variables. One may verify that evaluating 
(14.27)) on the relevant internal kinematic configuration reproduces eqs. ( 14 . 2 8 j) and (14.29)) . 

The component of the flying-squirrel cut of fig. 12(d) that we evaluated shows that this cut 
contributes to all two-loop scalar functions . The same is true for the other cuts with this 
topology but with cyclicly-permuted external legs. The sum over cyclic permutations may 
be reorganized in terms of permutations of a single "even" spin coefficient R\4& + -R413 which 
multiplies three different integrals of the type 1^ with different assignments of external legs. 
We will use this presentation in the following section. 

V. THE TWO-LOOP SIX-POINT SUPERAMPLITUDE 

We determined the four- dimensional cut-constructible even parts W^' ^ of the two- loop 
six-point NMHV superamplitude, 



W^ Mav = W\^ u= ^ + Wj , (5.1) 
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as explained in the previous section, by analyzing the cuts shown in fig.[2]or next-to-maximal 
versions of them. We obtained an explicit expression for the remaining part, W^' 1 *, cut- 
constructible only in D-dimensions, by comparing the results of D-dimensional and four- 
dimensional cut calculations. We have carried out the calculation without assuming a specific 
(possibly overcomplete) basis of two-loop integrals and found that the integrals listed in 
figs. H] and |8] are necessary and sufficient to saturate the cuts of the even part of the amplitude 
through O(e ). For comparison, and because of changes in the labeling of these i nteg rals 
with respect to the original calculations of the two-loop six-point MHV amplitude 19], we 
also present it in our labeling. 
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IS 




FIG. 8: /i-integrals entering the 2-loop 6-point amplitudes. The arrow on the external line indicates 
leg number 1. 



A. The NMHV amplitude 

The four- dimensional cut-constructible even part of all six-point two-loop amplitudes is 
built out of a sum of the 16 integrals shown in fig. H] 

S (2) ' D=4 (1) = \cj {1 \e) + c 2 I^\e) + \c^\e) + \cj^{e) + c 5 I^(e) + c^\e) 
+ \ (c 7a P- 2 /( 7 )( e ) + o^P- 1 JW(c) + c 7c I^(e)) + \c^\e) + c 9 lW(e) 

+ l -c l4 I^\t) + ic 15 /( 15 )(e) + c 16 / 16 (e) . (5.2) 

The coefficients q, which differ between the MHV and the NMHV amplitudes, are func- 
tions of external momenta and the numerical coefficients are symmetry factors reflecting the 
symmetries of each integral under cyclic permutations of external legs. 
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The functions W^' D 4 are constructed by summing S^ ,D 4 over the sets of permutations 



Si in eqs. (12. 9p . ( 12 . 10[) and (12. lip that map each superinvariant (R i+3 ^i +2 + -^+^+5) into 
itself: 



W, 



(2),D=4 



- y s®> 

8 ^ 



D=4, 



a) + 0(e) 



(5.3) 



Of the overall factor of 1/8, a factor of 1/4 emerges from the calculation of the unitarity 
cuts and a factor of 1/2 is due to our choice of normalization. 

The coefficients Cj in the identity permutation entering S^' D=4 in eq. (I5.2p are: 



Cl 



S 123 S34S 6 i + S123S234S345 



— S123S234S12S45 — S123S345S23S56 + 2Si2S23<S45<S56 

C3 = -5123(51235345 — S12S45) 

C5 = — S 12 Si23S234 
Cla — — ■5l23( s 345 s 234 — S6IS34) 
C7c = —S 123 (•5234-5 345 — -561 S34) 

C9 = Sl23 s 45"556 

c ll = — -556-561 -5123 
2 

c 13 = -5l23 s 61 
C15 = 

C17 = not necessary 



c 2 = 2Si 2 S 2 3 

2 

C4 — S123S34 

c 6 — S6I-S12S123 

Cjb — 2Si 2 3S34S6l 

c 8 = (5.4) 

c 10 — ■556-Sl23 s 345 

Cl2 = — s 123 (-5123-5345 ~~ S12.545) 

C14 = 2s| 4 Si23 

Ci6 = 2Si2-S34-5i23 

c 18 — |-5l23(2S34S61 — -5234-5345) 



In dimensional regularization, the six-point two-loop (and quite likely all higher-point 
higher-loop) amplitudes receive contributions from integrals — collectively referred to as 
"/x-integrals" - - whose integrand vanishes identically when evaluated in four dimensions. 
The integrals shown in fig. [8] are of this type, where /i p and [i q denote the (— 2e) components 
of the loop momenta. As noted in ref. 19], the integral J'- 17 -' vanishes identically as e — > 0; we 
will therefore ignore it in the following. To determine the contributions of such integrals we 
compare the result of the four- dimensional cut calculation with that of the .D- dimensional 
cuts and find that the even part of the amplitude also contains the terms, 



(2),M 



\aeSi J o-eSiUS 2 us 3 



(5.5) 



The coefficients Cj bear certain similarities to the corresponding coefficients in the MHV 
amplitude. 
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B. The MHV amplitude 



For completeness, and because of differences of notation from ref . 19[ , we also present the 
integrand of the even part of the MHV amplitude. The four-dimensional cut-constructible 
part is given by 

1 

16 



M (2),D=4 

where the coefficients Cj in the identity permutation are given by 



J2 S®> D = 4 (a) + 0(e) 

t7G<SiU5 2 U<S 3 



(5.6) 



C\ = 


Sl23 (•5l2'545' s 234 + S23 s 56 s 345 


c 2 = 


2S23 S 12 




+ S123(S34S61 - S234S345)) 






c 3 = 


s 123( s 345' s 123 — S45S12) 


C 4 = 


S 34'S?23 


c 5 = 


■Sl2( s 234 s 123 — 2S23S56) 


c 6 = 


— S61 s 12'Sl23 


C7a = 


S123(S234S345 — S34S61) 


C7b = 


-4s 34 S 6 iSi23 


Cjc — 


s 123( s 234'5345 — S34SQ1) 


c 8 = 


2512(5345^123 — S12S45) 


c 9 = 


•S45 s 56 s 123 


ClO — 


S56(2Si2S45 — S123S345) 


Cll = 


S61 s 56 s 123 


C12 = 


s 123( s 345 s 123 — s 12-545) 


Cl3 = 


-S123S6I 


C14 = 





Cl5 = 





Cl6 = 





Cl7 = 


— 2Si23<S345(s234<S345 — S61S34) 


CIS = 


2Si2(Si23S234S345 — Si 2 



(5.7) 



+25345(5123455234 + S23 s 56 s 345 

The //-integral contribution is 

r(2),M _ 1 



-2Si2(s23S56S345 + S34S61S123) 



Ma 
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E 



cr 



(5.8) 



o-e5iU5 2 uS3 

As mentioned previously, starts at (9(e) [19] and thus does not contribute through 

O(e ). 



C. A Comparison of the MHV and NMHV Amplitudes 

A direct inspection of the integrals in fig. H] and of their coefficients in eq. (15.41) re- 
veals that the even part of the two-loop six-point NMHV amplitude is a sum of pseudo- 
conformal integrals. This is similar to the MHV amplitude, for which the four- dimensional 
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cut-constructible part has a similar property 19], as may also be seen by directly inspecting 
the coefficients listed in eq. (|5.7p . 

This is perhaps not completely surprising in light of the argument presented in section UTTl 
that all four- dimensional cuts can be reproduced by cuts of pseudo-conformal integrals. This 
structure does not guarantee, however, that the even part of the amplitude is dual conformal 
invariant, even after infrared divergences are removed appropriately. We return to this point 
in the next section. 

The structure of the NMHV amplitude is quite similar to that of the MHV amplitude, 
with only subtle differences in the values of the coefficients. Two of the integrals that did 
not contribute to the MHV amplitude — J^ 14 ^ (e) and 1^ (e) — enter in the NMHV ampli- 
tude with nonvanishing coefficient; similarly, a topology that exists in both amplitudes - 
1^ (e) — appears in the NMHV amplitude with an additional pseudo-conformal numerator. 
Moreover, an integral that contributes to the MHV amplitude — 1^ (e) — disappears from 
the NMHV one. We note also that a perfectly valid integral — I^ 15 \e) — appears in neither 
the MHV nor NMHV amplitudes. It would be interesting to understand the significance of 
this observation. 

The properties of the MHV and NMHV amplitudes differ from those observed in the 
four-point amplitudes through five loops: 

• All pseudo-conformal integrals appear with relative weights of ±1 or [ol. [37!^ 3Q ] . 



An integral appears with coefficient zero if and only i: 



taking its external legs off shell and taking e — > 16]. 



the integral is unregulated after 



It has been proposed that the signs ±1 of the contributing integrals can be understood 



by the requirement of cancelling unphysical singularities 80]. 

It would undoubtedly be interesting to understand the generalization of these features to 
higher-multiplicity amplitudes. 

Although the amplitudes have very similar structures, the ratio of the NMHV six-point 
superamplitude to the MHV one does not appear to exhibit a transparent organization. This 
is due to the different structures of the permutation sums that contribute to the independent 
factors in the amplitudes. Indeed, each function W^' D ~ 4 and W^' D ~ 4 contains only a sum 
over the permutations in the set Si while the functions M^' D=4: and M^' D=4: contain sums 
over all twelve permutations Si U S 2 U S 3 . 
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VI. DUAL CONFORMAL INVARIANCE 



The explicit calculation in section IIVI of the two-loop NMHV six-point superamplitude 
shows that it indeed has the structure anticipated in section IIHI We obtained explicit 
integral representations for the scalar functions V4y , summarized in the previous section. 

As it is the case with all massless theories in four dimensions, the amplitude is infrared 
divergent; to examine the dual conformal properties of the amplitude it is necessary to 
isolate these divergences. Based on the universality of infrared singularities and their ex- 



ponentiation DHKS proposed 24J that these divergences be removed by simply dividing by 
the MHV amplitude. This ratio (I2.35P was conjectured to be dual conformally invariant. 
An alternative method was described in section lITTj see eq. (I3.16p . 

Since the functions Wf have a natural decomposition into four- dimensional and D- 
dimensional cut-constructible contributions, the functions Cf introduced in equation (13. 9p 
inherit a similar decomposition 

w (l) =w (l),D=4 + w (l),» _^ c (l) ^ c (l),D=4 + c p, fl l=1 ^ (6 J) 

At one loop, the /i-integral contribution cf^'^ vanishes in the limit e — > 0. Because of 
infrared divergences, they nevertheless give rise to nontrivial contributions at two loops in 
the ratio with the MHV amplitude. The functions Cf contain terms that vanish as e — > 0: 

Cf(e) = ^ (1) (6)-M«(6) 

Cf\e) = Wl 2 \e) - M™(e) - M« (e)(HA {1) (e) - M«(e)) . (6.2) 

As explained in sect. Ill C\ in the limit e — > 0, C[ 1 \e) reduces to dual conformal invariant 
functions very closely related to the V® defined in ref. [24j. The higher-order terms in e 
were not considered in ref. [24 1; we take them as determined by the amplitude through the 
relation flUD between W t (l) and Cf ] . 

The infrared- divergent terms in both (e) and (e) have the usual form 

6 

2c~ 



Div — 4D-^i)" £ • ( 6 - 3 ) 



Thus, C\ (e) is manifestly finite; moreover, the only divergent contribution arising from the 
last term in C„- (e) is due to the overall factor of M^(e). 



44 



A. Terms Requiring D-dimensional Cuts 



At one loop, the //-integrals, requiring consideration of I)- dimensional cuts, yield only 
terms of O(e) in both W i (e) and M^(e). This allows us to isolate the //-integral contri- 
bution in eq. (16.21) : 

Cf ^ = - wf^DivJ - + M^Div^ + 0(e) , (6.4) 

where we used the universality of the one-loop infrared divergences (16. 3p and kept only the 
terms that have nontrivial divergent and finite parts. For example, we dropped the terms 
in eq. (16. 2p coming from the finite part of the overall factor (e) in the last term. 

The last two terms in the equation above contain information already available in the 
MHV amplitude. Indeed, this exact combination appears in the iteration of the //-integrals 
for this amplitude 



M (2U = M^Div^ . (6.5) 



Thus, Cf^^ is given by, 



c (2u = _ w n# my w ; ( 6 . 6 ) 

with given by eq. (15.51) . 

This expression for Wf ^ may be further simplified by making use of the special prop- 
erties of the hexabox integral discussed in section IV. A of ref. [3], in particular eq. (4.6): 

^ U) ^ = -l(- 8 u)" 1_,/ta M • ( 6 - 7 ) 

Thus, ),M can be expressed exactly in terms of one-loop integrals, albeit in six dimensions. 
Moreover, using the fact that the one-loop hexagon integral is invariant under cyclic per- 
mutations of external legs, can be expressed in terms of the massless six-dimensional 
hexagon integral with a coefficient given by the universal divergent part of one-loop ampli- 
tudes: 



^Dh4 1} / he V] ^123(2534^61 - ^234*345) + °( € ) " ( 6 " 8 ) 
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The four terms in each sum are in fact equal. 

The //-integral contribution to the one-loop six-point NMHV amplitude is not yet avail- 
able in the literature. Information on its structure may be obtained by analyzing a two- 
particle cut of the //-integral contribution to the two-loop NMHV superamplitude. Dixon 
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and Schabinger [25| have evaluated such a cut directly; quite surprisingly, they find that it 
can be organized in terms of the same R invariants as the four-dimensional cut-constructible 
terms. The //-integrals' contribution to the even part of the one-loop six-point NMHV am- 
plitude is, 



^/ hC V] S 123(2S 3 4S 6 1 - S234S345) • (6.9) 
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Combining this with eqs. (16.81) and (I6.6P immediately shows that 

limd 2) ' M = . (6.10) 

e->0 

In other words, the complete /i-integral contribution to the six-point two-loop NMHV ampli- 
tude is completely accounted for by extracting an overall factor of the MHV superamplitude. 

Through similar manipulations it is possible to show that the remainder-like functions 
R & 1 introduced in eq. (13.161) do not receive any //-integral contributions. Indeed, directly 
expanding eq. (I3.16P to 0(a 2 ) we find that, 



l(W l {l \e)y + f 2 (e)W l W (2e) 



(6.11) 



Identifying the //-integral contributions to each of the terms on the right hand side and using 
the universality of infrared divergences implies that 

R<gj<» = W^(e) - Div^wf ] (e) = cf U + 0(e) . (6.12) 

It therefore follows from equations (I6.10p that R^^ does not receive any finite /i-integral 
contributions. 

The same, however, cannot be said about the /i-integral contribution to the odd part of 
the amplitude. Indeed, repeating the steps that lead to equations (16. 2p we find that the 
coefficients of the parity-odd quantities {Ri + ^,i,i+2 — Ri,i+z,i+b) are, 

Cf ) = Wf ] - Div^ Wf 1 . (6.13) 

' — (2) 

While the /i-integral contributions to W\ are given in terms of the hexabox integral or, 
equivalently in terms of the six-dimensional hexagon integral, their contributions to wf^ 
are given in terms of a restricted set of the one- mass pentagon integrals |25(] . This suggests 
that, for the odd part of the superamplitude, the /i integrals cannot be cleanly separated 
from the four-dimensional cut-constructible terms. 
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B. Numerical Evaluation of the Amplitude 



In order to further analyze the properties of the two-loop six-point NMHV amplitude, 
we turn to a numerical evaluation of the two-loop integrals. Thanks to the results described 
in the previous section, we may focus on the four- dimensional cut-constructible part of the 
amplitude. The task of evaluating the integrals is simplified substantially by the fact that 



all of them have already been evaluated at several distinct kinematic points in 



19|. We have 



evaluated additional kinematic points using the package MB [45| and the same Mellin-Barnes 
parametrization of integrals that was used in the calculation of the MHV amplitude. Apart 
from testing the symmetry properties of the amplitude, this calculation also verifies the 
expected universality of two-loop infrared divergences. Viewed differently, a successful test 
of the universality of the infrared divergences is a strong indication of the completeness of 
the cut construction described in previous sections. 

We choose Euclidean kinematics for all configurations of external momenta. As in the 
calculation of the MHV amplitude, the symmetries of the momentum configuration, 

K {0) : s i4+1 = -1 , s i>i+lii+2 = -2 , (6.14) 

make it particularly useful, as all cyclic permutations or external legs yield the same value 
for all integrals. This implies that all functions W^' D ~ 4: are equal for all i = 1,2, 3. Using 



the values of the integrals collected in the Appendix B of 



19| we find 



W d=*( K (o)\ = i + a ( + 5.27682 + 8.73314e + 8.11147e 2 ) 

e 2 



9 / 9 14.5967 25.3014 ±0.0043 \ „ 

+ a ( — j — + 21.064 ± 0.002 ) + O a 3 . 6.15 

V2e 4 e l e / 

Where they are not explicitly included, the errors do not affect the last quoted digit. We have 
used the error estimated reported by CUBA [81] . In general we found the errors to be reliable, 
giving an accurate measure of the number of trustworthy digits. In some contributions, 
however, we found them to be underestimated, invariably in the presence of small integrals 
with a fast- varying integrand. In such cases, when CUBA reports a large x 2 > we take the 
average value of the integrals to be the central value and quote the variation of the integral 
under changes of sampling points as the error estimate. The issue presumably involves 
integration regions missed because of special properties of the integrand. 

(2) 

As discussed previously, the construction of the functions C\ requires keeping higher 
orders in the small-e expansion of the one-loop amplitude. For the MHV amplitude, we use 
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(2) 

the expression in terms of the iterated one-loop amplitude and the remainder function R 6 . 
For the point we find 

Ci(a,e,K {Q) ) = l + a(0.783676 + 1.10087e + 0.07507e 2 ) (6.16) 
+ a'( - a003 ^ a ° 043 - (2.412 ±0.002) - /#>(*<»>)) + O(a'). 

We note that the residue of the simple pole in e vanishes within errors, as it should. We 
have confirmed this property for all the other kinematic points 9 . 

From equation (16.151) we can also find the value of the remainder-like functions R^.- 
introduced in equation (I3.16P at the point 

rW = _i.430 ±0.002 . (6.17) 

Similarly to the error quoted for Ci(a, e, K^), the error of R^j is completely inherited from 
that of Wf\ 

We have evaluated the amplitude at another kinematic point (denoted by K^> in \}\ ) 
related to by dual conformal transformations as well as two other points related to each 
other but unrelated to K^: 

: sia = -0.7236200, s 23 = -0.9213500, s 34 = -0.2723200, s 45 = -0.3582300, 
s 56 = -0.4235500, s 61 = -0.3218573, Si 23 = -2.1486192, s 234 = -0.7264904, 
s 345 = -0.4825841, 
K (3) : s iA+1 = -1, s 123 = -1/2, s 234 = -5/8, s 345 = -17/14, 

: s 12 = -2, s 23 = -4, s 3i = -2, s 45 = -14/17, s 56 = -4/5, s 61 = -56/85, 
Si,i+i,i+2 = -1, (6-18) 

In listing the kinematic points we attempted to preserve the notation for the points used in 
19]. We have collected our results for the values of Cf^ and R^j in tables fl] and [III The 



ref. 



three dual conformal ratios 

/ \ ( -Si2S 45 S 23 S 56 S 34 S 6 i \ 

(mi,m 2 ,m 3 )=( , , ) (6.19) 

v -5i 23 S 34 5 S 234 S 4 56 S 34 5S56l/ 

for the kinematic points are listed in the second column of these tables. 



9 We have also verified analytically the cancellation of infrared singular terms through 0(e 2 ). The complete 
cancellation of infrared-singular terms was shown analytically by G. Korchemsky (private communication). 
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TABLE I: Comparison of conformally-related kinematic points. C\ are the finite parts of the 

(2) 

ratios C, = W^/Mq at two-loops. i?g is the two- loop remainder function of the six-point MHV 
amplitude. 



kinematic pt. 


(ni,n 2 ,n 3 ) 


C? + Rf 


cf +4 2) 






U' 4' 


-2.413 ± 0.002 


-2.413 ±0.002 


-2.413 ±0.002 


K W 


U' 4' 4^ 


-2.359 ± 0.048 


-2.375 ±0.025 


-2.380 ± 0.033 




/28 16 112 \ 
U7' 5 ' 85 I 


14.426 ± 0.003 


12.614 ±0.004 


11.697 ±0.009 




(28 16 112 \ 
U7> 5 ' 85 > 


14.439 ± 0.078 


12.614 ±0.035 


11.727 ±0.145 



(2) 

TABLE II: Comparison of the remainder-like functions Rq.- at conformally-related kinematic 
points. 



kinematic pt. 


(ui,u 2 ,u 3 ) 


p(2) 


p(2) 
^6;2 


p(2) 




(1 1 I) 

\ 4 ) 4 > 4 7 


-1.431 ±0.002 


-1.431 ±0.002 


-1.431 ±0.002 




(I I I) 

\ 4 ' 4 ' 4 7 


-1.377 ±0.048 


-1.393 ±0.025 


-1.397 ±0.033 


K m 


(28 16 112 \ 
V 17' 5 ' 85 / 


5.413 ± 0.003 


4.749 ± 0.004 


4.602 ± 0.009 


K <6) 


(28 16 112 \ 
V 17 ' 5 ' 85 / 


5.427 ± 0.078 


4.749 ± 0.035 


4.633 ±0.145 



As mentioned previously, we left the MHV remainder function R\ unevaluated. Its dual 



12| 



conformal invariance 



(2), 



(6.20) 



implies that the equality within errors of the relevant entries of table [J extends to an equality 

(2\ 

of the functions C\ . Alternatively, we could have evaluated the remainder function from the 
analytic expression found in ref . 82] , the integral representation in ref . 83J , or the simplified 

(2) (2.) 

form in ref. 84]. The results we obtain thus suggest that C\ and R G .l are functions solely 
of the conformal cross-ratios, that is, that they are indeed invariant under dual conformal 
transformations. 
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VII. SUMMARY, CONCLUSIONS AND SOME OPEN QUESTIONS 



The maximally supersymmetric gauge theory in four dimensions is an ideal testing ground 
for probing the properties of gauge theories at both weak and strong coupling. The large 
degree of symmetry makes perturbative calculations tractable to relatively high orders while 
its string-theory dual provides powerful tools for understanding its strong-coupling behavior. 
Its hidden symmetries yield additional constraints that go beyond their initial connection 
to the integrability of the dilatation operator of the theory. 

In this paper we have computed the parity-even part of the two-loop six-point NMHV 
amplitude using generalized unitarity in superspace. We showed that the result is invari- 
ant under dual conformal transformations, after removal of universal infrared divergences 
(including terms arising from 0(e) contributions at one loop, computed by Dixon and Sch- 
abinger |25j). The dual conformal invariant content may be organized in several different 
ways. The exponentiation of both the infrared divergences and of the collinear splitting am- 
plitudes suggest the introduction of certain remainder-like functions which, similarly to the 
remainder function for MHV amplitudes, are functions only of the conformal cross ratios. 
We have shown that, to all orders in perturbation theory, it should be possible to reconstruct 
the remainder-like functions by evaluating certain triple-collinear splitting amplitudes. 

Several interesting issues related to the calculation described here, and to the structure 
of the perturbative expansion of the theory and its strong coupling expansion remain to be 
clarified. 

Through the AdS/CFT correspondence, Alday and Maldacena argued that, to leading 
order in their strong-coupling expansion, all planar scattering amplitudes with fixed number 
of external legs are essentially identical up to perhaps a rational function of momenta and 
polarization vectors. Our calculation and arguments show that the weak-coupling structure 
of the six point amplitude involves at least six different spin factors dressed with scalar 
and pseudo-scalar functions to all orders in the weak-coupling expansion. Reconciling this 
structure with the AdS/CFT considerations remains an important open problem, which 
appears to require that, in the strong-coupling limit, the remainder-like functions introduced 
in section II I II have a very simple relation to the MHV remainder function. 

Inspired by strong-coupling considerations, several groups showed that MHV amplitudes 
have a close relation to certain light-like polygonal Wilson loops order-by-order in weak- 
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conpfing perturbation theory, first in expfieit ca.cu,ation S MM, and very recent.y, ca 8 t 
in a more general setting [201 ] . A similar relation for non-MHV amplitudes remains an 
open question; our numerical results provide check points for future calculations in this 
direction. The Wilson-loop formulation of the six-point MHV amplitude led to the analytic 
evaluation 
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84| of the remainder function at two loops. It seems likely that a Wilson- 



loop formulation of NMHV amplitudes will allow a similar evaluation of the remainder-like 
functions characterizing this amplitude. 

A direct comparison of the integrands of the six-point MHV and NMHV amplitudes at 
two loops reveals that certain integrals appear in one but not the other, while the contribut- 
ing integrals enter with numerical coefficients that do not conform with the effective rules 
inferred from four-point amplitudes. Moreover, one perfectly valid pseudo-conformal inte- 
gral does not appear in either one of these amplitudes. It would be interesting to develop 
a better understanding of this pattern of numerical coefficients. Evaluation of higher-point 
NMHV amplitudes at two loops may help in this direction. 

In our calculation dual conformal invariance is obscured by the dimensional regulator; 
removal of infrared divergences is a crucial step in studying the dual conformal properties 
of scattering amplitudes. Using four and five-point amplitudes as testing ground, it was 
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that regulating the infrared divergences by a particular symmetry breaking of 
the gauge group makes dual conformal invariance more transparent. It would be interesting 
to repeat the calculation described in this paper as well as that of the two-loop six-point 
MHV amplitude in this framework. Apart from a better understanding of dual conformal 
invariance, such an endeavor would clarify the interpretation of /^-integrals as well as that of 
the parity-odd component of amplitudes. While we did not compute the parity-odd part of 
the six-point NMHV amplitude explicitly, its /x-integral contributions make it quite different 
from the parity-odd terms in MHV amplitudes with up to six external legs. 

A very interesting question relates to symmetries of scattering amplitudes. Up to anoma- 
lies introduced by the regulator, it is expected that dual conformal transformations leave 
scattering amplitudes invariant to all orders in perturbation theory. Again apart from 
anomalies due to the presence of a regulator, ordinary conformal invariance exhibits ad- 
ditional anomalies of a holomorphic type 86[, related to singular momentum configurations, 
already at tree level. An appropriate definition of a generating function for superamplitudes 
with variable number of external legs [87| allows this anomaly to be circumvented, and to 
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be realized on scatte ring amplitudes through the one-loop level. Korchemsky and Sokatchev 
have recently given [34J a general construction of conformal and dual conformal invariants 
(and hence of Yangian invariants). Other ap pro aches to the construction of Yangian invari- 



ants were discussed by Mason and Skinner 



and by Drummond and Ferro [89]. These 
by Arkani-Hamed et al. to represent 



invariants have expressions that were conjectured 
the leading singularities of scattering amplitudes to all orders in perturbation theory. All- 
order leading singularities have been derived directly by Bullimore, Mason, and Skinner (33}] . 
While the form of subleading singularities is not yet clear, this structure suggests that it 
may be possible to realize both symmetries (and their closure) at higher loops. The ex- 
tended algebra does not uniquely determine the S-matrix of M = 4 super- Yang-Mills theory, 



however 



34J. Unraveling its constraints on scattering amplitudes should prove fruitful. 



The study of non-MHV amplitudes at higher loops is only in its early stages. Our explicit 
calculation provides an example of such an amplitude. The structure of these amplitudes 
is substantially richer than that of MHV amplitudes. It seems likely that new and exciting 
properties as well as new calculational techniques are waiting to be discovered. 
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